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The  steady  aerodynamic  results  are  shown  in  the  form  of  upper  and  lower 
surface*  pressure  curves.  The  unsteady  aerodynamic  coefficients  are  ob¬ 
tained  for  various  values  of  low  reduced  frequencies  by  pitching  the  air¬ 
foils  about  the  quarter  chord  axis.  Unsteady  results  are  presented  as 
plots  of  unsteady  coefficients  versus  angle  of  attack  for  the  MBB  A-3 
airfoil  and  versus  Mach  number  for  the  CAST  7  and  TF-8A  airfoils.. 


At  design  Mach  number  of  0.765,  aerodynamic  and  corresponding  flutter  re¬ 
sults  are  obtained  for  the  MBB  A-3  airfoil  by  both  programs  at  various 
angles  of  attack.  The  angles  of  attack  considered  when  using  LTRAN2  are 
-0.2°,  0.0°,  0.2°,  0.42°,  0.6°,  and  0.8°,  respectively,  and  those  con¬ 
sidered  when  using  STRANS2/UTRANS2  are  -0.4°,  0.0°,  0.42°,  0.75°,  1.0°, 
1.2°,  and  1.3°,  respectively.  Results  are  presented  as  curves  of  flutter 
speed  and  corresponding  reduced  frequency  versus  angle  of  attack  for 
selected  values  of  aeroelastic  parameters.  The  effects  of  angle  of  attack 
on  flutter  speed  obtained  by  the  two  programs  are  compared  and  discussed. 

For  the  CAST  7  supercritical  airfoil,  aerodynamic  and  corresponding  flutter 
results  are  obtained  at  zero  mean  angle  of  attack  for  various  Mach  num¬ 
bers  by  using  LTRAN2.  The  Mach  numbers  considered  are  0.6,  0.625,  0.650, 
0.675,  0.70,  0.71,  and  0.72,  respectively.  Results  are  presented  as  curves, 
of  flutter  speed  and  the  corresponding  reduced  frequency  versus  Mach  num¬ 
ber  for  selected  values  of  aeroelastic  parameters.  The  transonic  dip 
phenomenon  is  observed  in  all  curves.  The  effect  of  each  aeroelastic 
parameter  on  the  curves  is  studied. 

For  the  TF-8A  wing  section  at  65.3%  semispan,  aerodynamic  and  corres¬ 
ponding  flutter  results  are  obtained  at  0.0°  and  -3.0°  angle  of  attack 
for  various  Mach  numbers  by  STRANS2/UTRANS2.  The  angle  of  attack  -3.0°  of 
the  present  wing  section  corresponds  to  zero  mean  angle  of  attack  of  the 
zero- twist  section  of  the  wing  in  the  wind  tunnel  experiment  conducted  by 
NASA.  The  Mach  numbers  considered  are  0.70,  0.74,  0.76,  0.77,  0.78,  0.79, 
and  0.80,  for  a  =  -3.0°,  and  0.70,  0.72,  0.74,  0.76,  and  0.78  for 
a  =  0.0°,  respectively.  Flutter  results  are  presented  in  the  form  of 
curves  of  flutter  speed  and  the  corresponding  reduced  frequency  versus 
Mach  number  for  selected  values  of  the  aeroelastic  parameters.  The 
transonic  dip  phenomenon  is  more  pronounced  in  the  curves  for  -3.0°  local 
angle  of  attack.  The  effect  of  each  aeroelastic  parameter  on  the  curves 
is  studied  and  discussed. 
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-  free  stream  velocity 

-  distance  between  mid-chord  and  pitching  axis  in  semi-chords 
positive  toward  the  trailing  edge 

-  S/mb,  distance  between  elastic  axis  and  center  of  mass  in 
semi -chords,  positive  toward  the  trailing  edge 

-  pitching  degree  of  freedom,  positive  in  nose-up  direction 

-  phase  angle  between  lift  force  and  plunging  displacement 

-  phase  angle  between  lift  force  and  pitching  rotation 

-  phase  angle  between  moment  and  plunging  displacement 

-  phase  angle  between  moment  and  pitching  rotation 

-  induced  angle  of  attack 

-  ratio  of  specific  heats 

-  h/c 

-  unsteady  perturbation  parameter 

-  flutter  eigenvalue 
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-  m/Trpb  ,  airfoil-air  mass  density  ratio 


-  h/b 

-  free  stream  air  density 

-  ratio  of  airfoil  thickness  to  chord  length 

-  disturbance  velocity  potential 

-  flutter  frequency  of  harmonic  oscillation 

-  (k^/m)  ,  uncoupled  plunging  frequency 

1  /z 

-  (k  /I  )  ,  uncoupled  pitching  frequency 

-  reference  frequency  set  equal  to  unity 
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SECTION  I 


INTRODUCTION 

Over  the  last  decade,  one  of  the  major  emphasis  of  the  research 
in  airfoil  technology  has  been  on  the  analysis  and  design  of  super¬ 
critical  airfoils.  Such  airfoils  have  proven  to  be  more  efficient 
than  conventional  airfoils,  particularly  in  the  transonic  regime.  In 
general,  supercritical  airfoils  have  higher  drag  rise  Mach  numbers  and 
higher  lift  coefficients  when  compared  to  equivalent  conventional  air¬ 
foils. 

Extensive  theoretical  and  experimental  studies  have  been  conducted 
to  prove  that  supercritical  airfoils  are  aerodynamically  more  efficient 
than  conventional  airfoils.  However,  comparisons  between  the  performance 
of  these  two  types  of  airfoils  based  on  aeroelastic  studies  have  begun 
only  recently.  The  purpose  of  this  report  is  to  study  the  flutter 
characteristics  of  supercritical  airfoils  in  the  transonic  regime. 

In  this  study,  three  supercritical  airfoils  were  considered:  (1) 

MBB  A-3  airfoil  designed  by  Messerschmi tt-Bolkow-Blohm  of  the  Federal  Re¬ 
public  of  Germany;  (2)  CAST  7  airfoil  designed  by  Dornier  GmbH  of  the 
Federal  Republic  of  Germany;  and  (3)  TF-8A  wing  section  at  the  65.3% 
semi  span  station  designed  by  NASA. 

1.  Development  of  Supercritical  Airfoils 

Early  literature  regarding  the  development  of  the  concept  of  super¬ 
critical  airfoils  can  be  found  in  References  1  and  2.  A  brief  report  on 


ci  recent  flutter  analysis  of  supercritical  airfoils  is  given  in 
Reference  3. 

In  Reference  4,  Whitcomb  discussed  research  conducted  on  super¬ 
critical  airfoils  at  NASA.  The  salient  aerodynamic  features  of  super¬ 
critical  airfoils  were  pointed  out.  Major  advantages  of  these  airfoils 
over  the  conventional  airfoils  were  illustrated.  Methods  of  incorporating 
supercritical  airfoils  in  swept  wings  were  discussed. 

In  Reference  5,  Whitcomb  discussed  flight  correlations  for  the  F-8 
supercritical  wing  configurations  developed  by  NASA.  Reasons  for  the 
differences  between  the  flight  and  wind  tunnel  results  for  pitching 
moment,  drag  polars,  drag  rise,  and  pressure  distribution  for  the  F-8 
wing  were  discussed.  The  primary  differences  were  attributed  to  Reynolds 
number  effects  and  wind  tunnel-wall  interference.  It  was  stated  that 
wall  interference  particularly  influenced  the  transonic  drag  rise. 

In  Reference  6,  Yu  presented  a  transonic  shock-free  wing  redesign 
procedure  by  using  a  fictitious  gas  method.  The  method  was  based  on 
the  full  potential  equation  and  local  redesign  of  the  wing  surface 
geometry  beneath  the  supersonic  region  to  produce  a  shock-free  flow.  Re¬ 
sults  for  a  redesigned  non-lifting  rectangular  wing  and  for  a  redesigned 
lifting  wing  of  ONERA  M6  planform  were  presented.  A  supercritical  wing 
configuration  obtained  by  this  method  showed  significant  drag  reduction 
and  improved  lift-drag  characteristics. 

A  detailed  description  of  the  design  and  analysis  of  supercritical 
wing  sections  was  given  by  Bauer,  Garabedian,  and  Korn  in  References  7,  8, 
and  9  with  Jameson  being  another  author  in  Reference  8.  In  Reference  7. 
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the  authors  presented  a  mathematical  theory  based  on  the  finite  difference 
scheme  for  the  design  of  supercritical  wing  sections.  They  presented  a 
computer  program  for  the  analysis  and  also  illustrated  the  method  by  ex¬ 
amples.  Several  shockless  airfoils  were  designed  and  tested  in  the  wind 
tunnel.  Their  theoretical  results  established  satisfactory  agreement  with 
the  experimental  results. 

In  Reference  8,  the  work  carried  out  in  Reference  7  was  modified  and 
was  presented  in  a  more  definitive  form.  The  authors  improved  the  design 
of  the  trailing  edge  by  which  a  lift  increase  of  15  to  20  percent  was  ob¬ 
tained.  The  improved  method  had  provision  to  handle  supersonic  as  well  as 
subsonic  free  stream  Mach  numbers,  and  to  capture  the  shock  wave  as  far 
back  on  an  airfoil  as  required.  Moreover,  the  method  led  to  an  effective 
three-dimensional  program  for  the  computation  of  transonic  flows  past  an 
oblique  wing.  Detailed  comparisons  were  made  with  experimental  data.  The 
comparisons  were  favorable. 

The  work  presented  in  Reference  9  was  a  sequel  to  the  two  earlier 
References  (7  and  8).  New  mathematical  techniques  for  the  design  and 
analysis  of  a  supercritical  airfoil  were  incorporated  into  the  computer 
code.  The  advanced  mathematical  approach  made  it  possible  to  assign  the 
pressure  as  a  function  of  the  arc  length  and  then  obtain  a  shockless  airfoil 
that  nearly  achieves  the  desired  distribution  of  pressure.  This  tool  en¬ 
abled  them  to  design  families  of  transonic  airfoils  more  easily  both  for 
airplane  wings  and  for  compressor  blades  in  cascades. 

The  iterative  process  employed  in  Reference  9  can  be  summarized  as 
follows:  (1)  The  airfoil  is  prescribed  ana  mapped  onto  the  unit  circle. 

The  free  stream  Mach  number  is  specified  and  either  the  coefficient  of  lift 
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or  the  angle  of  attack  can  be  prescribed;  (2)  The  flow  calculations  are 
executed  for  a  fixed  number  of  cycles;  (3)  A  boundary  layer  correction 
is  computed  and  added  to  the  original  airfoil  to  give  a  new  profile;  and 
(4)  The  profile  is  mapped  onto  the  unit  circle.  Steps  2  through  4  are 
repeated  till  a  satisfactory  profile  is  obtained.  The  authors  illustrated 
the  method  with  examples.  Results  compare  well  with  those  obtained  from 
experiment. 

2.  Wind  Tunnel  Studies  on  Supercritical  Airfoils 

Many  of  the  supercritical  airfoils  designed  have  been  tested  in  the 
wind  tunnel  for  aerodynamic  performance.  Some  of  these  wind  tunnel  tests 
are  discussed  in  this  section. 

In  Reference  10,  Bucciantini,  Oggiano,  and  Onorato  conducted  wind 
tunnel  tests  on  the  MBB  A-3  supercritical  airfoil.  Steady  state  surface 
pressure  distributions,  wake  and  boundary  layer  measurements  were  obtained 
in  two  different  wind  tunnels  by  using  the  same  model.  The  chord  length 
and  span  of  the  model  were  equal  to  0.127  meters  and  0.203  meters,  re¬ 
spectively.  Several  Mach  numbers  between  0.751  and  0.855  and  several 
angles  of  attack  between  1.12°  and  4.18°  were  considered.  Results  were 
presented  in  the  form  of  steady  pressure  distributions.  Measured  values 
of  the  lift  and  moment  coefficients  about  the  quarter  chord  axis  were 
also  presented.  These  results  can  serve  as  an  experimental  data  base  for 
comparison  with  computational  results  for  the  MBB  A-3  airfoil. 

In  Reference  11,  Stanewsky,  Puffert,  Muller,  and  Bateman  presented 
wind  tunnel  test  results  for  the  CAST  7  supercritical  airfoil.  Steady 
state  surface  pressure  distributions,  wake  and  boundary  layer  measurements 
were  obtained  in  three  independent  wind  tunnels.  Several  Mach  numbers 
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between  0.4  and  0.8  and  several  angles  of  attack  between  -2.0°  and  +5.0° 
were  considered.  Results  were  presented  in  the  form  of  steady  pressure 
curves.  Measured  values  of  drag,  lift  and  moment  coefficients, and 
boundary  layer  thickness  were  also  given.  These  results  can  serve  as  an 
experimental  data  base  for  comparison  with  computational  results  for 
the  CAST  7  airfoil. 

In  Reference  12,  Harris  and  Bartlett  conducted  wind  tunnel  experi¬ 
ments  on  a  NASA  supercritical  wing  research  airplane  model.  Experiments 
were  conducted  at  the  NASA-Langley  8-foot  transonic  pressure  tunnel  with 
Mach  numbers  ranging  from  0.25  to  1.00.  They  measured  the  steady  aero¬ 
dynamic  load  distributions  over  the  wing  and  rear  fuselage  of  the  mocel 
with  and  without  fuselage  area-rule  additions.  Pressure  measurements  ove-' 
the  surface  of  the  area-rule  additions  at  sideslip  angles  of  approximately 
-5°,  0°,  and  5°  were  also  included. 

In  Reference  13,  Davis  and  Malcolm  obtained  experimental  unsteady 
aerodynamic  results  for  a  conventional  and  a  supercritical  airfoil. 

The  two  airfoils  were  NACA  64A010,  a  10%  thick  airfoil  of  conventional 
shape,  and  NLR  7301,  a  16.5%  thick  supercritical  airfoil.  Results  were 
obtained  for  four  conditions:  (1)  conventional  airfoil  in  transonic  flow 
with  weak  normal  shock;  (2)  conventional  airfoil  in  transonic  flow  with 
stronger  shock;  (3)  supercritical  airfoil  at  experimentally  determined 
snock-free  design  condition;  and  (4)  supercritical  airfoil  at  off-design 
condition  with  stronger  shock.  Measured  chordwise  unsteady  pressure  time- 
histories  for  the  four  conditions  were  compared.  It  was  concluded  that 
although  an  oscillating  supercritical  airfoil  excites  more  harmonics,  the 
strength  of  the  airfoil  shock  is  the  more  important  parameter  governing 
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the  complexity  of  the  unsteady  flow.  Whether  they  are  conventional 
or  supercritical,  it  was  found  that  airfoils  which  support  weak  shocks 
induce  unsteady  loads  that  are  qualitatively  predictable  with  classical 
theories.  It  was  also  found  that  flows  with  strong  shocks  are  sensitive 
to  details  of  the  shock  and  boundary-layer  interaction  and  cannot  be  ade¬ 
quately  predicted.  Results  presented  can  provide  as  a  comparative  basis 
for  those  obtained  by  the  new  generation  of  Navier-Stokes-type  codes. 

3.  Developments  in  Unsteady  Transonic  Aerodynamics 

Due  to  the  fast  growth  of  digital  computers,  advances  in  the  numerical 
computational  methods  of  a  transonic  flow  field  around  oscillating  two- 
dimensional  airfoils  have  been  extensive.  Many  computer  programs  have  been 
developed  based  on  these  methods.  A  bibliography  of  these  developments  was 
given  by  Borland  in  Reference  14.  It  was  observed  in  References  15  and  3 
that  both  the  computer  programs  STRANS2  and  UTRANS2  developed  by  Traci, 
Albano,  and  Farr  (Reference  16)  and  the  LTRAN2  developed  by  Ballhaus  and 
Goorjian  (Reference  17)  can  be  efficiently  used  to  obtain  unsteady  transonic 
aerodynamic  data  for  conventional  and  supercriti cal  airfoils. 

The  computer  programs  STRANS2  and  UTRANS2  use  a  harmonic  analysis 
method  to  solve  the  two-dimensional ,  moderate-frequency,  small-disturbance 
transonic  equations.  The  unsteady  aerodynamic  equations  are  linearized 
with  respect  to  time.  Thus,  the  unsteady  solution  is  treated  as  a  small 
linear  perturbation  over  a  nonlinear  steady  state  solution.  The  position  of 
the  shock  is  fixed  at  the  steady  state  position  throughout  the  analysis. 

In  these  programs,  both  steady  and  unsteady  aerodynamic  equations  are  solved 
by  a  mixed  differencing  line  relaxation  procedure. 
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On  the  other  hand,  the  computer  program  LTRAN2  is  based  on  the  time 
integration  method.  This  program  solves  the  two-dimensional,  low-f requency , 
smal 1-di sturbance  transonic  equation.  The  solution  for  the  steady  part  is 
obtained  by  successive  line  over-relaxation  method  (SLOR).  The  solution  for 
the  unsteady  part  is  obtained  by  an  alternating-direction  implicit  finite 
difference  algorithm.  In  this  program, the  shock  is  allowed  to  move  during 
the  unsteady  analysis. 

In  Reference  18,  Houwink  and  Van  der  Vooren  modified  LTRAN2  by  ex¬ 
tending  the  program  for  high  frequency  computations.  High  frequency  terms 
were  added  to  boundary  conditions  and  pressure  computations.  However,  the 
same  basic  low  frequency  potential  equation  used  in  LTRAN2  was  again  con¬ 
sidered.  Modifications  showed  improvements  on  the  results  obtained  for  a 
flat  plate  at  M  =  0.7.  For  the  case  of  the  harmonic  pitching  oscillation 
of  a  NACA  64A006  airfoil  and  a  reduced  frequency  of  0.8  based  on  full  chord, 
they  found  significant  differences  relative  to  the  LTRAN2  results.  However, 
these  results  have  not  been  compared  with  other  transonic  codes. 

Recently,  Rizzetta  and  Yoshihara  (Reference  19)  developed  an  unsteady 
transonic  code  ExTRAN2  whi ch  does  not  have  any  low  frequency  assumptions. 

High  frequency  terms  were  retained  in  potential  equation,  boundary  condi¬ 
tions,  and  pressure  computations.  Viscous  effects  were  incorporated  in 
the  program  by  using  a  viscous  ramp  method.  Unsteady  pressure  and  co¬ 
efficients  were  computed  for  a  NACA  64A010  airfoil  at  M  ~  0.80.  It  was 
shown  that  these  results  compared  favorably  with  those  reported  in 
Reference  18.  Also  the  boundary  conditions  play  a  more  important  role 
than  does  the  differential  equation  in  obtaining  unsteady  aerodynamic  co¬ 
efficients.  In  addition,  viscous  interaction  was  found  to  have  the  expected 
important  effects. 
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4.  Developments  in  Transonic  Aeroelasticity  of  Supercritical  Airfoils 


Based  on  the  steady  aerodynamic  data,  the  supercritical  airfoils  have 
proven  to  be  more  efficient  than  conventional  airfoils  in  the  transonic 
regime.  However,  it  is  of  more  practical  interest  to  compare  the  unsteady 
aerodynamic  data  between  the  two  types  of  airfoils.  In  that  case,  aero- 
elastic  characteristics  of  the  two  types  of  airfoils  can  be  studied  and 
compared  and  more  useful  conclusions  can  be  drawn.  Studies  of  the  aero- 
elastic  characteristics  of  supercritical  airfoils  have  begun  recently. 

A  short  review  of  the  work  of  Farmer  and  Hanson  (Reference  20),  McGrew 
et  aj.  (Reference  21),  and  Ashley  (Reference  22)  was  given  by  the  authors 
in  Reference  23. 

In  Reference  3  and  23, the  present  authors  have  investigated  the  flutter 
characteri sties  of  a  MBB  A-3  supercritical  airfoil.  This  airfoil  is  one 
of  the  AGARD  standard  airfoils  suggested  for  aeroelastic  applications  of 
transonic  unsteady  aerodynamics  (Reference  24).  Results  presented  in 
References  3  and  23  are  briefly  discussed  here. 

The  steady  aerodynamic  results  were  obtained  by  both  5TRANS2  and 
LTRAN2  in  the  form  of  pressure  plots  on  upper  and  lower  surfaces.  The  un¬ 
steady  results  were  obtained  by  UTRANS2  (harmonic  method)  and  LTRAN2  (time 
integration  method)  in  the  form  of  lift  and  moment  aerodynamic  coefficients 
at  various  values  of  reduced  frequency. 

The  study  was  first  conducted  by  considering  a  case  at  design  Mach 
number  of  0.765  and  zero  mean  angle  of  attack.  Two  parallel  sets  of  results 
for  steady  pressure  distributions,  unsteady  aerodynamic  coefficients,  and 
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flutter  speed  curves  were  obtained  by  simultaneously  using  the  two  separate 
computer  programs.  The  two  sets  of  results  were  in  good  agreement. 

The  design  conditions  of  the  MBB  A- 3  supercri ti cal  airfoil  are: 

M  =  0.765,  angle  of  attack  u  =  1.3°,  and  steady  lift  coefficient  C  =  0.58.  From 

X 

the  steady  pressure  computations,  it  was  found  that  at  the  design  Mach  number 
of  0.765  the  angles  of  attack  required  to  produce  the  design  lift  coefficient 
of  0.58  were  0.75°  and  0.42°  by  STRANS2  and  LTRAN2,  respectively.  Hence, 
flutter  analyses  were  performed  separately  by  using  angles  of  attack  of  0.75° 
and  0.42°  for  the  two  respective  computer  programs.  Results  were  obtained 
and  discussed  in  a  fashion  similar  to  that  performed  in  the  case  of  zero  mean 
angle  of  attack.  Based  on  the  specific  values  assumed  for  the  aeroelastic 
parameters,  it  was  concluded  that  the  MBB  A-3  supercri ti cal  aeroelastic 
system  is,  in  general,  more  stable  at  the  equivalent  design  angles  of  attack 
than  at  zero  angle  of  attack. 

In  order  to  study  the  effect  of  camber  on  flutter  results,  a  MBB  A-3 
airfoil  without  camber  was  considered.  Flutter  analyses  were  performed  for 
design  Mach  number  of  0.765  and  at  angles  of  attack  of  0.75°  and  0.42°  for 
STRANS2/UTRANS2  and  LTRAN2,  respectively.  The  results  showed  that  the  camber 
has  a  beneficial  effect  in  increasing  the  flutter  speeds. 

A  parallel  set  of  results  was  also  obtained  for  a  NACA  64A010  conven¬ 
tional  airfoil  scaled  down  to  the  same  maximum  thickness-to-chord  ratio 
(8.97)  as  that  of  the  MBB  A-3  supercritical  airfoil.  The  aerodynamic  and 
flutter  results  were  compared  with  those  obtained  for  the  MBB  A-3  super¬ 
critical  airfoil,  with  and  without  camber.  It  was  found  that  the  scaled  down 
NACA  64 A0 1 0  airfoil  gives  aerodynamic  and  flutter  results  almost  identical 
to  those  for  the  MBB  A-3  airfoil  without  camber. 
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In  order  to  investigate  the  effect  of  Mach  number,  flutter  analysis  was 
also  carried  out  for  the  MBB  A-3  supercritical  airfoil  at  five  different  Mach 
numbers:  0.7,  0.72,  0.74,  0.765,  and  0.78.  The  angle  of  attack  was  assumed 
to  be  zero.  Flutter  curves  were  presented  as  plots  of  flutter  speed  and  the 
corresponding  reduced  frequency  versus  Mach  number  for  various  values  of 
different  aeroelastic  parameters.  The  "transonic  dip"  phenomenon  was  ob¬ 
served  and  the  comparisons  of  the  curves  by  the  two  computer  programs  were 
favorable. 

5.  Scope  of  the  Present  Study 

In  this  study,  flutter  analysis  was  conducted  for  three  supercritical 
airfoils:  (1)  MBB  A-3;  (2)  CAST  7;  and  (3)  TF-8A.  Aerodynamic  computations 
for  MBB  A-3  were  carried  out  by  both  the  STRANS2/UTRANS2  and  LTRAN2  transonic 
codes  simultaneously.  On  the  other  hand,  aerodynamic  computations  for  CAST  7 
and  TF-8A  were  carried  out  separately  by  LTRAN2  and  STRANS2/UTRANS2 ,  re¬ 
spectively.  For  all  the  airfoils,  two  degrees  of  freedom,  pitching  and 
plunging,  were  considered. 

Four  aeroelastic  parameters  were  considered  in  this  analysis:  (1)  airfoil- 
air  mass  density  ratio;  (2)  plunge-to-pitch  frequency  ratio;  (3)  position 
of  the  mass  center;  and  (4)  position  of  the  elastic  axis. 

a.  Flutter  Analysis  of  a  MBB  A-3  Supercritical  Airfoil 

This  study  is  a  sequel  to  the  study  reported  by  the  authors  in  Refer¬ 
ence  3.  Here,  the  effect  of  angle  of  attack  on  flutter  speed  was  studied. 
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The  steady  aerodynamic  results  were  obtained  by  both  STRANS2  and 
LTRAN2  in  the  form  of  pressure  plots  on  upper  and  lower  surfaces.  The 
unsteady  results  were  obtained  by  UTRANS2  and  LTRAN2  in  the  form  of  lift 
and  moment  coefficients  at  various  values  of  reduced  frequency. 

Because  of  the  different  limitations  of  the  two  programs,  separate 
values  of  angle  of  attack  were  considered  for  the  two  programs.  For  LTRAN2, 
six  angles  of  attack  between  -0.2°  and  0.8"  were  considered,  whereas  for 
STRANS2/UTRANS2 ,  seven  angles  of  attack  between  -0.4°  and  1.3°  were  con¬ 
sidered. 

Based  on  the  two  sets  of  aerodynamic  coefficients  obtained  for  various 
angles  of  attack,  flutter  analyses  of  the  MBB  A-3  airfoil  were  conducted. 

The  effect  of  angle  of  attack  on  flutter  speed  was  studied  for  various  values 
of  the  four  aeroelastic  parameters.  Results  obtained  by  the  two  programs 
are  compared  and  discussed. 

b.  Flutter  Analysis  of  a  CAST  7  Supercritical  Airfoil 

The  purpose  of  this  study  was  to  investigate  the  flutter  characteristics 
of  a  CAST  7  supercritical  airfoil.  This  configuration  is  one  of  the  AGARC 
standard  airfoils  suggested  for  aeroelastic  applications  of  transonic  un¬ 
steady  aerodynamics  (Reference  24). 

The  steady  and  unsteady  aerodynamic  results  were  obtained  by  LTRAN2. 
Seven  Mach  numbers  between  0.6  and  0.72  were  considered  at  zero  angle  of 
attack.  The  steady  results  are  presented  <n  the  form  of  plots  of  pressure 
curves  on  upper  and  lower  surfaces.  The  unsteady  results  are  presented  as 
plots  of  the  magnitudes  of  the  unsteady  coefficients  and  the  corresponding 
phase  angles  versus  Mach  number. 
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Based  on  the  unsteady  aerodynamic  coefficients  obtained  for  various 
Mach  numbers,  flutter  analysis  of  the  CAST  7  airfoil  was  conducted.  The 
effect  of  Mach  number  on  flutter  speed  was  studied  for  various  values  of 
the  four  aeroelastic  parameters.  Curves  of  flutter  speed  versus  Mach  number 
were  plotted.  The  transonic  dip  phenomenon  similar  to  that  obtained  for 
the  NACA  64A006  airfoil  (Reference  25)  and  the  MBB  A- 3  airfoil  (Reference 
3)  was  observed. 

c.  Flutter  Analysis  of  a  TF-8A  Wing  Section  at  the  65.3%  Semispan  Station 

The  purpose  of  this  study  was  to  investigate  flutter  characteristics 
of  a  TF-8A  wing  section  at  the  65.3%  semispan  station.  This  wing  was  de¬ 
signed  by  NASA  and  its  configuration  was  given  in  Reference  26.  Based  on 
experimental  data,  transonic  flutter  studies  on  this  wing  were  carried  out 
in  References  20,  21,  and  22. 

During  wind  tunnel  tests  on  the  TF-8A  wing  (Reference  26)  it  was  noticed 
that  because  of  twist,  the  angle  of  attack  at  the  65.3%  semispan  station 
was  -3.0°.  This  corresponded  to  a  zero  angle  of  attack  of  the  zero-twist 
section  of  the  wing.  Based  on  the  wind  tunnel  data  given  in  Reference  26, 
the  transonic  dip  phenomenon  was  found  for  the  TF-8A  wing  section  in  Re¬ 
ference  22.  Hence,  the  present  flutter  analysis  for  the  TF-8A  wing  section 
was  performed  at  a  mean  angle  of  attack  equal  to  -3.0°  as  well  as  zero 
degrees . 

The  steady  and  unsteady  aerodynamic  data  were  obtained  by  STRANS2  and 
UTRANS2,  respectively.  Seven  Mach  numbers  between  0.70  and  0.80  were  con¬ 
sidered  for  the  airfoil  at  -3.0°  mean  angle  of  attack  and  five  Mach  numbers 
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between  0.70  and  0.78  were  considered  for  the  airfoil  at  zero  mean  angle 
of  attack.  Steady  state  results  are  presented  in  the  form  of  pressure 
curves  on  the  upper  and  lower  surfaces.  Unsteady  results  are  presented  as 
plots  of  magnitudes  of  unsteady  coefficients  and  corresponding  phase  angles 
versus  Mach  numbers. 

Based  on  the  unsteady  aerodynamic  coefficients  obtained  for  these  Mach 
numbers  at  the  two  angles  of  attack,  flutter  analysis  of  the  TF-8A  wing 
section  was  conducted.  The  effect  of  Mach  number  on  flutter  speed  was 
studied  for  various  values  of  the  four  aeroelastic  parameters.  Curves  of 
flutter  speed  and  corresponding  reduced  freguency  versus  Mach  number  were 
plotted.  The  transonic  dip  phenomenon  was  observed. 


SECTION  II 


TRANSONIC  FLOW  EQUATIONS  AND  TWO  COMPUTATIONAL  METHODS 


Developments  of  numerical  procedures  for  obtaining  practical  solutions 
for  the  unsteady  flow  around  two-dimensional  airfoils  have  been  extensive. 

A  number  of  solution  procedures  have  led  to  successful  transonic  codes.  The 
computer  codes  based  on  time  integration  method  and  harmonic  method  are  two 
among  them.  In  this  section  the  aerodynamic  equations  employed  in  the  two 
codes  are  discussed. 

1.  Unsteady  Two-Dimensional  Flow  Equation  for  Transonic  Flow 

The  simplified  basic  aerodynamic  equation,  following  the  assumptions 
that  the  flow  is  two-dimensional,  inviscid,  transonic  (M^  =1),  and  that  the 
velocity  disturbances  are  small  compared  to  the  free  stream  velocity  U,  can 
be  deduced  from  the  general  equation  of  continuity  of  gas  dynamics  as 

k2  M2  <*..  +  2k  M2d>  .  =  V  <j>  +  <t  (1) 

c  «  vtt  c  »*xt  c  vxx  vyy  v  ' 

where  kc  =  Uc/U  is  the  reduced  frequency;  Mm  is  the  free  stream  Mach  number; 
<f>  is  the  disturbance  velocity  potential;  Vc  =  1  -  M2  -  (y  +  1  )m£  4>x;  m  is  a 
function  of  M  ;  and  y  is  the  ratio  of  specific  heats. 

CO 

In  deriving  the  above  equation,  the  coordinate  system  is  fixed  with 
respect  to  the  airfoil,  and  x  is  aligned  with  the  free  stream  direction.  The 
flow  is  defined  as  locally  subsonic  or  supersonic,  relative  to  the  fixed 
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coordinate  system,  for  Vc>0  or  Vc<0,  respecti vely.  A  measure  of  the  degree 
of  unsteadiness  is  given  by  the  reduced  frequency  kc  when  the  airfoil  is 
oscillating  periodically  with  a  frequency  w. 

Several  numerical  approaches  have  been  developed  for  the  solution  of 
transonic  flow  fields  governed  by  Equation  1. 

2.  Harmonic  Method 

In  the  harmonic  method  it  is  assumed  that  the  flow  field  for  some 
sinusoidal  body  motion  of  frequency  w  can  be  expressed  in  the  form 

4>(x,y,t)  =  <J)Q(x ,y)  +  e4>i  (x,y)e"lwt  +  c2f2(x,y)e1wt 

+  higher-order  terms  (2) 

where  $  is  the  disturbance  velocity  potential,  and  £  is  related  to  the 

amplitude  of  body  motion.  For  purely  subsonic  or  supersonic  flows,  the 

sinusoidal  motion  produces  a  sinusoidal  response  <j>  at  the  same  frequency 

and  all  higher-order  terms  are  zero.  This  is  not  true  in  the  transonic 

case,  in  which  a  higher  harmonic  content  in  4>  results  because  the  governing 

equations  are  nonlinear.  However,  if  the  amplitudes  of  motion  are  assumed 

2 

to  be  very  small  (e  <<  1),  terms  of  order  e  or  higher  can  still  be  neglected. 
Hence,  Equation  2  can  be  simplified  as 

♦  (x.y.t)  =  4>0(x,y)  +  e^U.yJe10’1  (3) 

In  Equation  3,  the  solution  for  depends  on  the  mean  steady-state 
solution  if.  .  In  this  equation,  it  is  implrr  that  the  unsteady  solution  is  a 
small  linear  perturbation  about  some  nonlinear  steady  state  solution.  In 
other  words,  the  unsteady  solution  is  linearized  with  respect  to  time. 


15 


Equation  3  has  the  advantage  that  4>-j  can  be  computed  using  essentially 
the  same  well-known  finite-difference  relaxation  algorithms  used  to  compute 
the  mean  steady  state  solution  <p  .  Based  on  this  approach,  Traci  et  a_l_. 
(Reference  16)  developed  computer  programs  STRANS2  and  UTRANS2  that  can  solve 
the  steady  and  unsteady  aerodynamic  equations,  respectively.  Because  of  the 
elliptic/hyperbolic  characteristic  of  the  governing  equations  for  the  steady 
and  unsteady  perturbation  velocity  potentials,  the  mixed  differencing  line 
relaxation  procedure  of  Murman  and  Cole  was  used  for  steady  and  unsteady 
computations.  These  programs  are  used  in  the  present  report  to  compute 
aerodynamic  data  for  flutter  analysis. 

3.  Time  Integration  Method 


In  the  harmonic  method,  a  time  linearized  assumption  is  used  to  obtain 
the  unsteady  solution.  Such  method  constrains  the  shock  wave  to  its  steady- 
state  position.  As  an  alternative  and  more  complete  solution,  Ballhaus  and 
Goorjian  (Reference  17)  proposed  to  use  the  time  integration  method  to  obtain 
an  unsteady  transonic  solution.  In  this  method,  the  shock  movements  can 
accurately  be  treated.  Also,  there  is  no  need  to  linearize  the  aerodynamic 
equations  with  respect  to  time  in  the  unsteady  part. 

As  a  further  simplification  of  Equation  1,  the  frequency  of  the  transonic 

2  2/3 

flow  can  be  assumed  as  low  so  that  kc  =  1  -  M^  -  t  <<  1 .  Equation  1  may 
then  be  reduced  to 


2k  Ml  =  V  4  +  <t 

c  ^xt  cvxx  \yy 

where  t  is  the  thickness-chord  ratio  of  the  airfoil. 


(4) 
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Equation  4  is  suitable  for  the  time  integration  approach.  This  approach 


is  based  on  the  finite-difference  scheme  that  integrates  Equation  2  in  time 
for  harmonic  aerodynamic  motions  until  the  transient  states  in  the  solution 
disappear  and  the  forces  become  periodic. 

Several  numerical  procedures  are  available  to  solve  Equation  4.  Among 
them  the  procedure  developed  by  Ballhaus  and  Goorjian  (Reference  17)  based  on 
the  alternate-direction  implicit  algorithm  has  been  proven  to  be  computat ional ly 
efficient  and  is  being  widely  used.  This  procedure  uses  a  conservative,  im¬ 
plicit  finite-difference  scheme  to  time-accurately  integrate  the  nonlinear, 
low-frequency,  transonic  small-disturbance  equation  as  defined  in  Equation  4. 

A  computer  code  LTRAN2  was  developed  based  on  this  procedure.  This  cede  can 
be  used  to  find  the  flow  field  solutions  for  airfoils  with  arbitrary  combina¬ 
tions  of  pitch,  plunge,  and  flap  deflections. 

In  order  to  comply  with  the  low-frequency  approximation  of  LTRAN2,  the 
unsteady  aerodynamic  coefficients  were  computed  for  a  reduced  frequency  range 
kc  <  0.2. 

Although  UTRANS2  retains  the  4>tt-term,  it  experiences  numerical  instability 
difficulty  at  higher  reduced  frequencies.  Hence,  in  this  study  the  reduced 
frequencies  considered  for  UTRANS2  were  also  limited  to  the  range  kc  <  0.2. 


4.  Procedures  for  Obtaining  Aerodynamic  Data 

A  description  of  the  procedures  for  obtaining  steady  and  unsteady  aero¬ 
dynamic  data  by  using  both  LTRAN2  and  STRANS2/UTRANS2  was  given  in  detail  in 
Section  IV  of  Reference  3.  The  reader  is  thus  referred  to  that  reference. 

In  that  section,  the  sizes  and  patterns  of  the  finite  difference  meshes 
at  various  stages  of  computation,  distributions  of  gridpoints  near  airfoils, 
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size  and  number  of  time  steps  required  for  convergence  in  integration  method, 
relaxation  parameters  and  convergence  tolerances  in  relaxation  metnod,  and 
relevance  of  all  these  parameters  to  computer  capacity,  computing  time,  and 
accuracy  were  discussed  in  detail. 
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SECTION  III 


AEROELASTIC  AND  TRANSFORMATION  EQUATIONS 


In  the  present  flutter  analysis,  the  airfoil  is  assumed  to  oscillate 
about  the  elastic  axis  with  two  degrees  of  freedom,  plunge  and  pitch.  The 
necessary  aeroelastic  equations  of  motion  and  the  transformation  equations 
for  aerodynamic  coefficients  are  discussed  in  this  section. 

1.  Aeroelastic  Equations  of  Motion 


The  parameters  and  sign  conventions  for  a  typical  airfoil  oscillating 
with  pitching  and  plunging  degrees  of  freedom  are  defined  in  Figure  1.  The 
following  assumptions  are  made  in  deriving  the  equations: 

a.  The  displacement  h  and  rotation  a  are  measured  with  respect  to 
the  mean  position  defined  by  the  steady  state  conditions. 

b.  The  airfoil  is  rigid. 

c.  The  amplitudes  of  oscillation  are  small. 

d.  The  principle  of  superposition  for  aerodynamic  forces  is  valid 
even  in  the  presence  of  shocks.  Discussion  and  justification  of 
this  assumption  were  given  in  References  27  and  28. 

Assuming  harmonic  oscillations  with  flutter  frequency  w,  the  final 
eigenvalue  equations  for  flutter  analysis  are 


[y^[M]  -  [A]] 


(b) 
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Figure  1.  Definition  of  Parameters  for  Two-D.O.F.  Aeroelastic  System 


-n 

where  f q  and  aQ  are  the  nondimensional  amplitudes  in  plunge  and  pitch  oscilla- 

2 

tions,  respectively;  p  =  m/npb  ,  the  airfoil-air  mass  density  ratio; 

=  u,b/U,  the  reduced  frequency. 

The  matrices  [M],  [A],  and  [K]  are  defined  as 


The  eigenvalue  A  is  a  complex  number  defined  as 

A  =  y(l  +  ig)u,2b2/U2  (7) 

where  g  =  g^  =  g^  is  the  structural  damping  coefficient  which  is  assumed 
to  be  small  and  of  the  same  order  for  both  plunging  and  pitching  modes. 

The  flutter  solution  is  obtained  when  g  is  found  to  be  zero.  In  the  transonic 
flutter  analysis,  to  compute  the  aerodynamic  coefficients  for  Equation  6b 
comprises  the  essential  task. 


f  r.in‘,f  oniidf.  iori  of  Aerodynamic  Coefficients 


V 


In  some  situations, the  aerodynamic  coefficients  originally  obtained 
with  reference  to  a  pitching  axis  need  to  be  transformed  to  those  with  re¬ 
ference  to  the  elastic  axis  for  aeroelastic  analysis.  The  transformation 
relations  can  be  derived  by  using  the  principle  of  superposition  of  airloads. 

Let  it  be  assumed  that  the  aerodynamic  coefficients  C‘  ,  C!  ,  C'  ,  and 

£6  la  m<5 

C1  for  point  O'  (see  Figure  1)  due  to  pitching  about  O'  have  been  obtained. 

It  is  desired  to  transform  these  coefficients  to  be  C  C„  ,  C  .,  and  C  for 

SL&  SLa  m6  ma 

another  pitching  axis  (elastic  axis)  0. 

Assuming  rigid  airfoil  and  small  amplitude  of  oscillations,  the  trans¬ 
formation  relationships  may  be  written  as 


where  s  =  (a^  -  x  )/2.  The  three  terms  in  the  parentheses  are  relatively 
small  when  compared  with  the  other  terms.  If  they  are  neglected.  Equations  8 
are  in  the  form  as  those  given  previously  by  Traci  et  al.  in  Reference  16. 


SECTION  IV 


FLUTTER  ANALYSIS  OF  A  MBB  A-3  SUPERCRITICAL  AIRFOIL 


This  study  is  a  sequel  to  the  investigation  reported  by  the  authors 

in  Reference  3.  The  effect  of  angle  of  attack  on  flutter  speed  of  the  MBB  A-3 

supercritical  airfoil  was  studied  at  design  Mach  number  0.765. 

Aerodynamic  data  were  obtained  by  both  LTRAN2  and  STRANS2/UTRANS2. 

The  airfoil  was  assumed  to  oscillate  with  two  degrees  of  freedom,  plunge 
and  pitch,  about  the  elastic  axis. 

Because  of  the  different  limitations  of  the  two  programs,  separate 
values  of  angle  of  attack  were  considered,  respectively.  For  LTRAN2,  six 

angles  of  attack  between  -0.2°  and  0.8°  were  considered.  For  STRANS2/UTRANS2 , 

seven  angles  of  attack  between  -0.4°  and  1.3°  were  considered. 

Based  on  the  two  sets  of  aerodynamic  data  obtained,  flutter  analysis  of 
the  MBB  A-3  airfoil  was  conducted.  The  effect  of  angle  of  attack  on  flutter 
speed  was  studied  for  selected  values  of  four  different  aeroelastic  parameters. 
Results  obtained  by  the  two  programs  are  compared  and  discussed. 

Throughout  this  study,  the  values  for  radius  of  gyration  r  and 

a 

reference  frequency  w  were  assumed  as  0.5  and  1.0,  respectively. 

1.  Airfoil  Configuration 


In  order  to  perform  accurate  aerodynamic  computations,  it  is  necessary 
to  use  accurate  data  for  the  airfoil  configurations,  especially  in  the  region  of 


the  leading  edge.  The  equations  for  one-segment  fitting  of  eight  AGARD/SMP 
standard  airfoil  configurations  were  provided  by  Dr.  J.  J.  Olsen  of  the 
Flight  Dynamics  Laboratory  (Reference  24). 

Zu(x)  =  c ( x )  +  0.5  t(x ) 

Zg,( x )  =  c(x)  "  °-5 

234  (9) 

t  ( x )  =  aQVnr  +  a1  +  a2x  +  a3x  +  a4x  +  a5x 

c ( x )  =  ^1  +  ^2*  +  ^3X  +  ^4X  +  ^5X 

where  Zu(x)  and  Z^(x)  are  the  equations  for  the  upper  and  lower  surfaces, 
respectively;  x  is  the  nondimensional  axis  with  leading  edge  at  x  =  0  and 
trailing  edge  at  x  =  1 ;  c(x)  is  the  camber;  and  t(x)  is  the  thickness. 
Reference  24  also  gives  two-segment  fitting  functions  which  would  be 
necessary  for  more  precise  description  of  the  geometry. 

The  coefficients  for  the  thickness  function  for  the  MBB  A-3  airfoil  are: 
a  =  0.2457807;  a,  =  0;  a~  =  -0.2470393;  a,  =  0.5556936;  a.  =  -1.1377743;  and 
ac  =  0.5833393.  The  coefficients  for  the  camber  function  are:  b,  =  0; 

0  I 

b2  =  0.1294408;  b3  =  -0.4206819;  b4  =  0.5516741;  and  b5  =  -0.2604330.  A  plot 
of  this  configuration  is  given  in  Figure  2. 

The  maximum  thickness-to-chord  ratio  of  this  airfoil  is  8.9%.  It  has 
a  blunter  nose  than  equivalent  conventional  airfoils  (Reference  3). 

2.  Results  Based  on  LTRAN2 

Aerodynamic  data  were  obtained  by  the  time  integration  program  LTRAN2 
for  the  MBB  A-3  airfoil  at  design  Mach  number  0.765.  The  values  for  mean 


MBB  A-3 


angle  of  attack  considered  were  -0.2°,  0.0°,  0.2°,  0.42';,  0.6°,  and  0.8°. 
respectively,  where  0.42°  corresponds  to  the  equivalent  design  angle  of 
attack  that  yields  the  design  lift  coefficient  of  0.58  (Reference  3). 

a.  Steady  Pressure  Curves 

Figure  3  shows  the  steady  pressure  distributions  on  the  upper  and 
lower  surfaces  for  the  MBB  A- 3  airfoil  for  seven  mean  angles  of  attack  ob¬ 
tained  by  the  successive  line  over  relaxation  method  of  LTRAN2. 

The  pressure  curves  of  the  upper  surface  are  seen  to  move  upward  with 
the  increase  in  angle  of  attack  whereas  the  pressure  curves  of  the  lower 
surface  move  downward  with  the  increase  in  angle  of  attack.  Shocks 
are  seen  to  exist  for  all  angles  of  attack.  The  shock  moves  toward  the 
trailing  edge  with  increasing  strength  as  the  angle  of  attack  is  increased. 

From  the  nature  of  the  changes  in  the  steady  pressure  curves  it  can  be 
expected  that  both  lift  and  moment  coefficients  increase  as  the  angle  of 
attack  increases. 

h.  Unsteady  Aerodynamic  Coefficients 

Four  unsteady  aerodynamic  coefficients,  C„t,  C.  ,  C  .,  and  C  ,  were 

to  mo  nia 

obtained  by  pitching  the  airfoil  about  the  quarter  chord  axis  at  M  =  0.765. 
The  angles  of  attack  considered  were  -0.2°,  0.0°,  0.2°,  0.42°,  0."°,  and 
0.8°.  For  LTRAN2,  0.8"  was  practically  the  highest  angle  of  attack  that 
could  be  considered  for  this  airfoil.  Reduced  frequencies  considered  were 
0.05,  0.10,  0.15,  and  0.20,  respectively.  The  results  are  given  in  Table  1. 


STEADY  PRESSURE  COEFFICIENT  (c„) 


Table  1.  Aerodynamic  Coefficients  for  MBB  A-3  Airfoil  for 
Various  Angles  of  Attack  at  M  =  0.765  by  LTRAN2 


Reduced  Frequency  (kr) 


of 

i  Attack 

0 

.05 

0 

.10 

0. 

15 

0 

20  : 

Real 

Imaq. 

Real 

Imaq. 

Real 

Imaq . 

Real 

Imaq. 

•  -0.2° 

0.160 

0.598 

0.429 

1 .036 

0.694 

1.362 

0.977 

1 .595 j 

|  0.0° 

0.171 

0. 636 

0.422 

1.100 

0.724 

1 .420 

0.881 

1  .728 

0.2° 

0.238 

0.671 

0.595 

1.095 

0.900 

1  .386 

1.199 

1.563) 

0.42° 

0.324 

0.728 

0.734 

1.130 

1  .058 

1.379 

1.355 

1.505- 

« 

0.6° 

0.452 

0.783 

0.917 

1.133 

1.247 

1.314 

1  .493 

1 .41 7  | 

0.8°  |  0.709 

0.830 

1.165 

1.105 

1 .441 

1.231 

1 .619 

1.311' 

-0.2°  |  11.965 

-3.206 

10.356 

-4.290 

9.078 

-4.625 

7.974 

-  H 

; 

-4.886; 

0.0° 

12.724 

-3.409 

11.001 

-4.223 

9.469 

-4.825 

8.642 

-4.403: 

0.2° 

13.419 

-4.752 

10.952 

-5.946 

9.237 

-5.999 

7.815 

-5.997; 

0 . 4  2 

14.555 

-6.480 

11.303 

-7.340 

9.195 

-7.055 

7.524 

-6.776' 

* 

0.6° 

15.662 

-9.043 

11.329 

-9.174 

8.763 

-8.316 

7.086 

-7.467; 

0.8° 

16.604 

-14.18 

11.053 

-11.65 

8.205 

-9.607 

6.554 

-8 . 093' 

i 

:  -0.2° 

-0.005 

-0.026 

-0.031 

-0.049 

-0.022 

-0.073 

-0.030 

-0.100 

1 

|  0.0^ 

-0.008 

-0.038 

-0.015 

-0.071 

-0.022 

-0.103 

-0.029 

-0.136! 

;  o.2° 

-0  021 

-0.054 

-0.052 

-0.089 

-0.077 

-0.118 

-0.107 

-0. 1 40) 

;  0.42° 

-0.043 

-0.085 

-0.092 

-0.127 

-0.137 

-0.152 

-0.180 

-0.162 

0.6° 

-0.078 

-0.120 

-0.158 

-0.158 

-0.215 

-0.174 

-0.268 

-0.1 64i 

0.8° 

-0.171 

-0.171 

-0.280 

-0.182 

-0.346 

-0.176 

-0.385 

-0.171: 

:  -0.2° 

-0.511 

0.109 

-0.48° 

0.131 

-0.489 

0.145 

-0.502 

0.149- 

i 

;  O.lT 

-0.763 

0.162 

-0.706 

0.150 

-0.686 

0.146 

-0.680 

0.  145 

0.2° 

-1.084 

0.416 

-0.892 

0.515 

-0.790 

0.513 

-0.699 

0 . 537 ! 

*  0.42r- 

-1.704 

0.868 

-1.272 

0.924 

-1 .014 

0.913 

-0.811 

0.901 ' 

0.6° 

-2.391 

1.553 

-1.579 

1.579 

-1 .158 

1  .430 

-0.820 

1.338i 

!  0.8° 

-3.409 

3.409 

-1  .820 

2.803 

-1.175 

2 . 306 

-0.857 
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The  unsteady  coefficients  C„x,  C,  ,  C  .,  and  C  are  also  plotted 
against  angle  of  attack  in  Figures  4,  5,  6,  and  7,  respecti vely .  The 
corresponding  phase  angles  are  also  shown  in  the  four  respective  figures. 

In  Figure  4,  the  lift  coeff i cient , C£6! due  to  plunging  increases 
gradually  with  increase  in  angle  of  attack  whereas  the  corresponding  phase 
angle  9^  decreases  with  the  increase  in  angle  of  attack.  With  the  increase 
in  reduced  frequency, |C  |  increases  whereas  the  corresponding  e-j  decreases. 

In  Figure  5,  the  lift  coefficient  | j  due  to  pitching  about  the 
quarter  chord  axis  and  its  corresponding  phase  angle  both  increase  with 
the  increase  in  angle  of  attack.  With  the  increase  in  reduced  frequency, 

1C  I  decreases  whereas  I  0oi  increases. 

In  Figure  6,  the  magnitude  of  the  moment  coefficient  |C  6|  about  the 
quarter  chord  axis  due  to  plunging  increases  with  the  increase  in  angle  of 
attack,  whereas  the  corresponding  phase  angle  decreases  with  the  in¬ 
crease  in  angle  of  attack.  With  the  increase  in  reduced  frequency,  |Cmfij 
increases  and  6^  decreases. 

In  Figure  7,  the  magnitude  of  the  moment  coefficient  |C  |  about  the 
quarter  chord  axis  due  to  pitching  increases  with  the  increase  in  angle  of 
attack  whereas  the  corresponding  phase  angle  decreases  with  increase  in 
angle  of  attack.  With  the  increase  in  reduced  frequency,  |C  |  increases 
and  decreases. 

c.  Flutter  Results 

Based  on  the  unsteady  aerodynamic  coefficients  obtained  in  Table  1  for 
six  angles  of  attack,  flutter  analysis  was  performed.  The  emphasis  was  to  in¬ 
vestigate  the  effect  of  angle  of  attack  on  flutter  speed  for  various  values  of 
aeroelastic  parameters. 
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Figure  5.  Effect  of  Angle  of  Attack  on  Lift  Coefficient  due  to 
Pitching  for  MBB  A- 3  Airfoil  by  LTRAN2 . 
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ICIENT  (PITCH) 


Figure  8  shows  the  curves  for  flutter  speed  and  the  corresponding  re¬ 
duced  frequency  versus  the  mean  angle  of  attack  for  five  different  values  of 
the  position  of  mass  center  (x  =  0.25,  0.3,  0.4,  0.5,  and  0.6,  respectively). 
The  values  for  the  airfoil-air  mass  density  ratio  y,  plunge-to-pi tch  frequency 
ratio  uj  /u)  ,  and  the  position  of  the  elastic  axis  a,,  were  assumed  as  100, 

0.1,  and  -0.5,  respectively . 

It  is  seen  in  Figure  8  that  the  flutter  speed  increases  as  the  mass 

center  approaches  the  elastic  axis  or  as  x  becomes  smaller.  For  smaller 

a 

values  of  x^  and  higher  values  of  a,  difficulty  is  encountered  in  obtaining 
a  flutter  solution.  For  a  constant  value  of  x  ,the  flutter  speed  increases 

a 

with  the  increase  in  angle  of  attack.  Such  increasing  trends  are  more  rapid 

at  higher  angles  of  attack  and  at  lower  values  of  x  .  These  changes  are  the 

direct  effect  of  changes  in  the  unsteady  coefficients  shown  in  Table  1. 

Figure  9  shows  the  curves  for  flutter  speed  and  the  corresponding  reduced 

frequency  versus  angle  of  attack  for  three  different  values  of  airfoil-air 

mass  density  ratio  (y  =  100,  200,  and  300,  respectively).  The  values  for  the 

other  aeroelastic  parameters  x  ,  a,  ,  and  w./w  were  assumed  as  0.5,  -0.5, 

an  n  a 

and  0.1 ,  respectively. 

In  Figure  9,  the  flutter  speed  increases  with  the  increase  in  airfoil- 
air  mass  density  ratio.  Reduced  frequency  decreases  with  the  increase  in 
y.  The  flutter  speed  increases  with  the  increase  in  angle  of  attack.  The 
rate  of  increase  is  higher  at  higher  angles  of  attack. 

Figure  10  shows  the  curves  for  flutter  speed  and  the  corresponding  re¬ 
duced  frequency  versus  angle  of  attack  for  two  different  values  of  plunge- 
to-pitch  frequency  ratio  (w^/u^  =  0.1  and  0.2,  respectively).  The  values 
for  the  other  aeroelastic  parameters  y,  x  ,  and  a^  were  assumed  as  100,  0.5, 
and  -0.5,  respecti vely . 
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Figure  8.  Effect  of  Angle  of  Attack  on  Flutter  Speed  for  Five  Positions 
of  Mass  Center  for  MBB  A-3  Airfoil  by  LTRAN2. 
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REDUCED  FREQ.  (k„)  FLUTTER  SPEED 


ANGLE  OF  ATTACK  (a) 


Figure  10.  Effect  of  Angle  of  Attack  on  Flutter  Speed  for  Two  Values 
of  Plunge-to-Pitch  Frequency  Ratio  for  MBB  A-3  Airfoil  by 
LTRAN2. 
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It  is  seen  in  Figure  10  that  the  flutter  speed  curve  for  the  case 

u v / u  =  0.1  is  higher  than  that  for  the  case  of  wu/w  =  0.2.  The  reduced 
ha  ha 

frequency  curve  for  the  former  case  is  lower  than  that  for  the  latter. 

The  flutter  speed  curve  for  oju/uj  =  0.1  increases  with  the  increase  in 

angle  of  attack.  The  increase  is  more  rapid  at  higher  angles  of  attack. 

On  the  other  hand,  the  flutter  speed  curve  for  wu/w  =  0.2  first  increases 

ha 

lightly  with  the  increase  of  angle  of  attack  (a  <  0.0°)  and  then  decreases 
with  the  increase  in  angle  of  attack. 

Figure  11  shows  the  curves  for  flutter  speed  and  the  corresponding  re¬ 
duced  frequency  versus  the  angle  of  attack  for  three  different  positions  of 
elastic  axis  (a^  =-0.1,  -0.3,  and  -0.5,  respectively).  The  mass  center  was 

fixed  at  mid-chord  and  values  for  y  and  w./u  were  assumed  as  100  and  0.1, 

n  a 

respectively . 

It  is  seen  in  Figure  11  that  when  the  position  of  elastic  axis  is  varied 
from  ah  =  -0.1  to  a^  =  -0.5,  there  is  not  much  change  in  flutter  speed. 

All  the  curves  corresponding  to  =  0.1  in  Figures  8  to  11  show  that 

flutter  speed  increases  with  the  increase  in  angle  of  attack.  Such  increase 
becomes  more  rapid  at  higher  angles  of  attack. 

3.  Results  Based  on  STRANS2/UTRANS2 

In  this  section,  the  aerodynamic  data  were  obtained  for  the  MBB  A-3 
supercritical  airfoil  by  the  harmonic  analysis  programs  STRANS2/UTRANS2  at 
design  Mach  number  M  =  0.765.  The  values  for  angle  of  attack  considered 
were  -0.4°,  0.0°,  0.42°,  0.75°,  1.0°,  1.2°,  and  1.3°,  respecti vely .  The 
value  a  =  0.75°  corresponds  to  the  equivalent  design  angle  of  attack  that 
yields  a  design  lift  coefficient  of  0.58  (Reference  3). 


REDUCED  FREQ.  (kc)  FLUTTER  SPEED 


a.  Steady  Pressure  Curves 


Figure  12  shows  the  upper  and  lower  surface  steady  pressure  curves  for 
the  MBB  A- 3  airfoil  for  eight  mean  angles  of  attack  ranging  from  -1.2°  to 
1.3°  at  M  =  0.765  obtained  by  STRANS2. 

In  Figure  12,  the  pressure  curves  for  the  upper  surface  move  upward  with 
the  increase  in  angle  of  attack  whereas  the  pressure  curves  for  the  lower 
surface  move  downward  with  the  increase  in  angle  of  attack.  A  shock  appears 
at  an  angle  of  attack  of  about  0.0°  and  grows  stronger  with  the  increase  in 
angle  of  attack.  Similar  results  were  obtained  by  LTRAN2  in  Figure  3.  However, 
the  shock  in  Figure  12  remains  almost  in  the  same  position  (x/c  =  0.5) 
for  all  angles  of  attack.  On  the  other  hand,  the  shock  in  Figure  3  moves 
towards  the  trailing  edge  with  the  increase  in  angle  of  attack.  For  the 
same  flow  conditions,  the  shock  obtained  by  LTRAN2  is,  in  general,  stronger 
than  that  obtained  by  STRANS2. 

From  the  nature  of  the  changes  in  the  steady  pressure  curves  shown  in 
Figure  12,  it  can  be  expected  that  both  lift  and  moment  coefficients  increase 
with  increase  in  angle  of  attack. 

b.  Unsteady  Aerodynamic  Coefficients 

Based  on  the  steady  state  pressure  curves  shown  in  Figure  12,  four  un¬ 
steady  aerodynamic  coef f icients ,  C..,  C.  ,  C  .,  and  C  .were  obtained  by 

£o  Zol  mo  m  a' 

UTRANS2  for  the  MBB  A- 3  supercritical  airfoil.  The  airfoil  was  assumed  to 
pitch  about  the  quarter  chord  at  design  Mach  number  M  =  0.765.  Seven  angles 
of  attack  (-0.4",  0.0°,  0.42°,  0.75°,  1.0°,  1.2",  and  1.3")  were  considered. 
Coefficients  were  computed  at  five  reduced  frequencies  (0.0,  0.05,  0.10.  0.15, 
and  0.20). 
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Figure  12.  Effect  of  Angle  of  Attack  on  Steady  Pressure  Curves  for 
MBB  A- 3  Airfoil  by  STRANS2 . 


Table  2  shows  the  real  and  imaginary  parts  of  the  four  unsteady  aero¬ 


dynamic  coefficients.  The  four  coefficients  C„r,  C„  ,  C  and  C  and  their 

1  £6  £a  m<5  ma 

corresponding  phase  angles  e ^ ,  e2 ,  o^,  and  e^,  are  plotted  against  the  angle 
of  attack  in  Figures  13,  14,  15,  and  16,  respectively. 

In  Figure  13,  the  magnitude  of  the  lift  coefficient  due  to  plunge 
C,  ,  does  not  appear  to  vary  much  with  the  change  in  angle  of  attack. 

Similar  behavior  is  observed  for  the  corresponding  phase  angle  9^.  With  the 
increase  in  reduced  frequency,  increases  whereas  decreases. 

In  Figure  14,  the  magnitude  of  the  lift  coefficient  due  to  pitch  |C?6| 
about  the  quarter  chord  axis  does  not  appear  to  vary  much  with  the  change  in 
angle  of  attack.  Similar  behavior  is  observed  for  the  corresponding  phase 
angle  e^.  With  the  increase  in  reduced  frequency,  |  C£  J  decreases  whereas 
i  Op [  increases. 

In  Figure  15,  the  magnitude  of  the  moment  coefficient  about  the  quarter 

cnord  1C  !  due  to  plunge  increases  gradually  with  the  increase  in  angle  of 
1  ma' 

attack.  The  corresponding  phase  angle  0^  gradually  decreases  with  the  increase 
in  angle  of  attack.  The  values  of  | C  j  are  higher  for  higher  reduced  fre¬ 
quencies,  whereas,  the  curves  for  do  not  show  a  definite  trend. 

In  Figure  16,  the  magnitude  of  the  moment  coefficient  about  the  quarter 
chord  jC  i  due  to  pitch  about  the  quarter  chord  increases  with  the  increase 
in  angle  of  attack  a.  The  corresponding  phase  angle  decreases  with  the 
increase  in  angle  of  attack.  With  the  increase  in  reduced  frequency,  j C J 
increases  for  «  less  than  about  0.5"  and  decreases  for  a  greater  than  about 
0.5",  whereas  increases  for  all  a-values  considered. 

The  magnitudes  of  the  moment  coefficients  seem  to  reach  a  local  maximum 
at  <  ~  1.2°  whereas  the  corresponding  phase  angles  seem  to  reach  a  local 
minimum  at  this  angle  of  attack. 
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Table  2.  Aerodynamic  Coefficients  for  MBB  A-3  Airfoil  for  Various  Angles  of  Attack 
at  M  =  0.765  by  UTRANS2 


Reduced  Frequency  (kc) 


Real 

Imag. 

0.104 

0.492 

0.0! 

0.106 

C.  499 

0.0 

o.m 

0.508 

0.0 

0.119 

0.557 

0.0 

0.094 

0.513 

0.0 

0.089 

0.521 

0.0 

0.091 

0.512 

1.3°  110.376  0.0  9.345 

1.2°  *10. 590  0.0  9.514 

1.0°  1 1 0.709  0.0  9.610 

.  7 5°j  1 1 . 388  0.0  10.171 

. 42° 11.168  0.0  9.992 

0.0°  11.221  0.0  10.141 

-.4°  11.129  0.0  9.979 


Real  Imag. 
”07273  0.835 
0.281  0.844 
0.285  0.854 
0.279  0.906 
0.286  0.881 
0.260  0.900 
0.281  0.881 

8.254  -2.451 
8.389  -2.521 
8.463  -2.567 
8.932  -2.818 
8.757  -2.687 
8.843  -2.801 
8.780  -2.605 


0.492 
0.501 
0.508 
0.539 
0.536 
0.523 
0.521 

7.260  -2.718 
7.360  -2.808 
7.419  -2.853 
7.757  -3.144 
7.677  -2.948 
7.667  -3.060 
7.745  -2.857 


.087 
.100  0.7 
.172  0.7 
.138  0.7' 


-0.007  -0.044  -0.016  -0.077  |-0  031  -0.106  -0.051 

-0.006  -0.043  -0.016  -0.076  -0.029  -0.106  -0.050 

-0.005  -0.040  -0.012  -0.071  -0.023  |  -0.100  -0.040 

-0-003  -0.035  -0.003  -0. 060  {-0.009  j  -0.090  -0.020 

0.0  -0.024  0.002  -0. 047  0.002  -0.072  0.0 

0.003  -0.015  0.010  -0.031  0.017  -0.053  0.020 

0.003  -O.OI2!  0.010  -0.026  0.018  -0.044  0.025 


0.901  0.0 

0.899  0.0 

0.824  0.0 

0.691  0.0 

3.510  0.0 

0.311  0.0 

0.248  0.0 


0. 

079 

1 

O 

.715 

0. 

.073 

-0. 

.715 

0. 

,040 

-0. 

.672 

0. 

.023 

-0. 

.597 

1 

0 

.087 

-0. 

.476 

0 

I 

,148 

-0. 

,340 
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Figure  14.  Effect  of  Angle  of  Attack  on  Lift  Coefficient  due  to 
Pitching  for  MBB  A- 3  Airfoil  by  UTRANS2 . 
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Figure  15.  Effect  of  Angle  of  Attack  on  Moment  Coefficient  due  to 
Plunging  for  MBB  A-3  Airfoil  by  UTRANS2. 
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Figure  16.  Effect  of  Angle  of  Attack  on  Moment  Coefficient  due  to 
Pitching  for  MBB  A-3  Airfoil  by  UTRANS2. 


When  comparing  the  results  obtained  by  LTRAN2  in  Figures  4  to  7  and 
those  obtained  by  UTRANS2  in  Figures  13  to  16,  it  is  difficult  to  draw  con¬ 
clusions  consistent  among  all  figures.  However,  in  both  cases  the  moment 
coefficients,  j  Cffl6j  and  |CmJ  ,  increase  with  the  increase  in  angle  of  attack. 
Also,  all  the  coefficients  in  both  cases  vary  with  the  reduced  frequency  in 
a  similar  way  except  for  a  >  0.5°  in  Figure  16. 

c.  Flutter  Results 

Based  on  the  unsteady  aerodynamic  coefficients  shown  in  Table  2  for 
seven  angles  of  attack,  a  flutter  analysis  of  the  MBB  A-3  supercritical  air¬ 
foil  was  performed.  The  emphasis  was  to  investigate  the  effect  of  angle  of 
attack  on  flutter  speed  for  various  values  of  aeroelastic  parameters. 

Figure  17  shows  the  curves  for  flutter  speed  and  the  corresponding  re¬ 
duced  frequency  versus  angle  of  attack  for  six  values  of  the  position  of  mass 
center  (x^  =  0.1,  0.2,  0.3,  0.4,  0.5,  and  0.6,  respectively).  The  values 
for  the  airfoil-air  mass  density  ratio  y,  plunge-to-pitch  frequency  ratio 
wh/wa>  and  the  position  of  elastic  axis  a^,  were  assumed  as  100,  0.1,  and 
-0.5,  respectively . 

In  Figure  17,  the  flutter  speed  increases  rapidly  as  the  mass  center 
approaches  the  elastic  axis  or  as  xa  becomes  smaller.  For  x^  =  0.1,  diffi¬ 
culty  is  encountered  in  obtaining  a  flutter  speed.  The  reduced  frequency  in¬ 
creases  with  the  increase  in  x  .  For  a  constant  value  of  x  .the  flutter 

a  a 

speed  increases  with  the  increase  in  angle  of  attack.  Such  increasing  rates 

are  more  rapid  in  the  neighborhood  of  a  =  0.6°  to  1.1°  and  at  lower  values  of 

x  . 
a 

Figure  18  shows  curves  for  the  flutter  speed  and  the  corresponding  re¬ 
duced  frequency  versus  angle  of  attack  for  three  different  values  of 
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Figure  17.  Effect  of  Angle  of  Attack  on  Flutter  Speed  for  Six 
Positions  of  Mass  Center  for  MBB  A-3  Airfoil  by 
STRANS2/UTRANS2. 
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Figure  18. 
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Effect  of  Angle  of  Attack  on  Flutter  Speed  for  Three 
Values  of  Airfoil-Air  Mass  Density  Ratio  for  MBB  A-3 
Airfoil  by  STRANS2/UTRANS2. 
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airfoil-air  mass  density  ratio  (p  =  TOO,  200,  and  300).  The  values  for 

the  other  aeroelastic  parameters  x  ,  aK,  and  w,/w  were  assumed  as  0.5, 

anna 

-0.5,  and  0.1,  respectively. 

In  Figure  18,  the  flutter  speed  increases  with  the  increase  in 
airfoil-air  mass  density  ratio.  The  reduced  frequency  decreases  with  the 
increase  in  p.  The  flutter  speed  increases  with  the  increase  in  angle  of 
attack.  The  increase  rate  appears  to  be  higher  near  higher  angles  of  attack. 
This  is  due  to  the  rapid  changes  in  the  moment  coefficients  at  higher  angles 
of  attack  (see  Table  2).  However,  the  flutter  speed  curves  appear  to 
peak  near  a  -  1.2°. 

Figure  19  shows  the  curves  for  flutter  speed  and  the  corresponding  re¬ 
duced  frequency  versus  angle  of  attack  for  two  different  values  of  plunge- 
to-pitch  frequency  ratio  =  0.1  and  0.2,  respectively).  The  other 

aeroelastic  parameters  p,  x^,  and  ah  were  assumed  as  100,  0.5,  and  -0.5, 
respectively . 

In  Figure  19,  the  flutter  speed  curve  for  =0.1  is  higher  than 

that  for  co.  / d  =  0.2.  The  corresponding  curves  for  reduced  frequency  show 
n  a 

an  opposite  trend.  Both  flutter  curves  increase  with  the  increase  in  angle 
of  attack. 

Figure  20  shows  the  curves  for  the  flutter  speed  and  the  corresponding 
reduced  frequency  versus  angle  of  attack  for  two  different  values  of  position 
of  elastic  axis  (a^  =  -0.1  and  -0.5,  respectively).  The  mass  center  was  fixed 
at  mid-chord  and  the  values  for  p  and  <^/ioa  were  assumed  as  100  and  0.1, 
respectively. 
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REDUCED  FREQ.  (kc)  FLUTTER  SPEED  (U/ba>a) 


Figure  19.  Effect  of  Angle  of  Attack  on  Flutter  Speed  for  Two 

Values  of  Plunge- to-Pi tch  Frequency  Ratio  for  MBB  A- 3 
Airfoil  by  STRANS2/UTRANS2. 


RED.  FREQ.  (kc)  FLUTTER  SPEED  (U/bwa) 


In  Figure  20,  the  flutter  curve  for  a^  =  -0.1  is  slightly  higher  than 
that  for  ah  =  -0.5.  The  flutter  speed  increases  with  the  increase  in  angle 
of  attack.  The  rate  of  increase  is  higher  at  higher  angles  of  attack. 

Such  behavior  may  be  due  to  the  rapid  changes  in  the  moment  coefficients 
at  higher  angles  of  attack. 

Figures  17  to  20  show  that  the  flutter  speed  increases  with  the  in¬ 
crease  in  angle  of  attack.  Comparing  Figures  8  to  11  with  Figures  17  to  20, 

it  is  seen  that  for  the  case  of  u  /u  =  0.1,  all  the  curves  for  flutter  speed 

n  a 

increase  with  the  increase  in  angle  of  attack.  However,  for  wu/w  =  0.2,  the 

n  a 

flutter  speed  trend  with  angle  of  attack  is  different  for  the  two  sets  of 
aerodynamics  (LTRAN2  and  UTRANS2). 
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SECTION  V 


FLUTTER  ANALYSIS  OF  A  CAST  7  SUPERCRITICAL  AIRFOIL 

The  CAST  7  supercritical  airfoil  was  designed  by  Dornier  GmbH  of  the 
Federal  Republic  of  Germany  and  is  one  of  the  AGARD  standard  airfoils  sug¬ 
gested  for  aeroelastic  applications  of  transonic  unsteady  aerodynamics 
(Reference  24).  In  this  study,  flutter  characteristics  of  the  airfoil  in 
the  transonic  regime  are  investigated  by  using  LTRAN2. 

The  airfoil  was  assumed  to  oscillate  with  two  degrees  of  freedom,  plunge 
and  pitch,  about  the  elastic  axis.  The  aerodynamic  data  were  obtained  by 
pitching  the  airfoil  about  the  quarter  chord  axis  with  zero  mean  angle  of 
attack.  Seven  Mach  numbers  between  0.6  and  0.72  were  considered. 

Based  on  the  aerodynamic  data,  a  flutter  analysis  of  the  CAST  7  airfoil 
was  conducted  and  the  effect  of  Mach  number  on  flutter  speed  for  various 
values  of  aeroelastic  parameters  was  studied. 

Throughout  this  study,  the  values  for  the  radius  of  gyration  r ^  and  the 
reference  frequency  were  assumed  as  0.5  and  1.0,  respectively. 

1.  Airfoil  Configuration 

The  configuration  of  the  CAST  7  airfoil  was  obtained  by  Equation  9. 

The  coefficients  for  one-segment  fitting  of  the  thickness  (provided  by 
Dr.  J.  J.  Olsen)  are:  aQ  =  0.3315970;  a]  =  0.0;  a?  =  -0.2393648;  a3  =  0.3331678; 
a.  =  -1.0215039;  and  ac  =  0.6030865.  The  coefficients  for  the  camber  function 
are:  b]  =  0.0;  b?  =  0.1216041;  b3  =  -0.5012554;  b4  =  0.8749108;  and  b5  = 
-0.4983959.  A  plot  of  the  configuration  is  given  in  Figure  21. 
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It  is  seen  in  Figure  21  that  the  maximum  thickness-to-chord  ratio  is 
11. 9k- .  The  airfoil  has  a  fairly  blunt  nose  when  compared  to  conventional 
airfoils. 

2.  Steady  Pressure  Curves 

CAST  7  is  a  fairly  thick  supercritical  airfoil  (11.92)  when  compared 
to  other  supercritical  airfoils  such  as  MBB  A-3  (8.92).  LTRAN2  is  based 
on  inviscid  and  small  disturbance  transonic  theory  which  is  valid  for  thin 
airfoils  (i.e.  <  122).  Hence,  it  was  of  interest  to  compare  present  LTRAN2 
steady  pressure  results  with  those  obtained  by  experiment. 

Figure  22  shows  two  sets  of  steady  pressure  curves,  one  obtained  by 
LTRAN2  and  the  other  by  experiment  in  Reference  11.  The  values  for  Mach 
number  and  mean  angle  of  attack  were  equal  to  0.602°and  0.5°,  respectively. 

The  two  sets  of  curves  appear  to  have  similar  overall  shapes.  LTRAN2 
gives  higher  pressure  coefficients,  thus  higher  lift.  The  discrepancy 
appears  to  be  larger  on  the  lower  surface  near  the  trailing  edge. 

Reasons  for  the  discrepancy  between  the  two  sets  of  curves  may  be 
attributed  to:  (1)  viscous  effects  which  are  not  considered  in  LTRAN2; 

(2)  small  disturbance  theory  used  in  LTRAN2;  and  (3)  possible  errors  intro¬ 
duced  due  to  both  numerical  and  experimental  difficulties. 

Figure  23  shows  the  steady  pressure  curves  for  the  upper  and  lower 
surfaces  for  eight  different  Mach  numbers  (0.60,  0.625,  0.65,  0.675,  0.70, 
0.71,  0.72,  and  0.73,  respecti vely)  obtained  by  LTRAN2.  With  increase  in 
Mach  number,  the  upper  surface  pressurt  curves  change  faster  than  the  lower 
surface  pressure  curves.  A  shock  develops  at  FI  '  0.7.  When  Mach  number  in¬ 
creases,  the  shock  moves  toward  the  trailing  edge  with  increasing  strength. 


STEADY  PRESSURE  COEFFICIENT  (c 


CAST  7  1  M ■ 0.602  ,  ot  =0.5° 
- LTRAN2,  - EXPERIMENT 


0  0.2  0.4  x/c  0.6  0.8  1.0 


Figure  22.  Comparison  between  Steady  Pressure  Curves  by  LTRAN2 
and  Experiment  for  CAST  7  Airfoil. 


58 


STEADY  PRESSURE  COEFFICIENT  (c 
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Figure  23.  Effect  of  Mach  Number  on  Steady  Pressure  Curves  for  CAST  7 
Airfoil  by  LTRAN2 . 
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(•'.  1  h*>  Itn.h  number  -i[i|<r ouches  0.73,  the  shock  approaches  the  trailing  edge. 
It  appears  that  Mach  number  0.7?  is  practically  the  highest  Mach  number  that 
can  be  considered  in  the  present  analysis. 


3.  Unsteady  Aerodynamic  Coefficients 

Four  unsteady  aerodynamic  coefficients,  C^,  ,  Cffl(5 ,  and  Cma,  were 

obtained  for  the  CAST  7  airfoil  by  pitching  the  airfoil  about  the  quarter 
chord  axis  with  zero  mean  angle  of  attack.  Seven  Mach  numbers,  0.60,  0.625, 
0.650,  0.675,  0.70,  0.71,  and  0.72,  were  considered.  Reduced  frequencies 
assumed  were  0.05,  0.10,  0.15,  and  0.20,  respectively.  The  coefficients 
are  presented  in  Table  3  and  plotted  in  Figures  24  to  27  against  Mach  num¬ 
ber.  Also  plotted  are  the  corresponding  phase  angles. 

In  Figure  24,  the  magnitude  of  the  lift  coefficient  due  to  plunge  jC^j 
gradually  increases  with  increase  in  Mach  number,  whereas,  the  corresponding 
phase  angle  e-j  decreases.  With  the  increase  in  reduced  frequency,  ]C  e6j  in¬ 
creases  and  G-|  decreases.  The  changes  in  both  jC£(S)  and  are  more  rapid 
above  M  =  0.70. 

In  Figure  25,  both  the  magnitude  of  the  lift  coefficient  [C£6|  due  to 
pitch  about  the  quarter  chord  axis  and  the  corresponding  phase  angle  | e ^  | 
increase  as  Mach  number  increases.  When  the  reduced  frequency  increases, 
jC£6|  decreases, whereas, | |  increases.  All  curves  take  a  sharper  increase 
as  M  becomes  larger  than  approximately  0.7. 

In  Figure  26, the  magnitude  of  the  moment  coefficient  |C  |  about  the 
quarter  chord  axis  due  to  plunge  increases  with  increase  in  Mach  number. 

The  corresponding  phase  angle  first  increases  up  to  M  <  0.68  and  then 
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PHASE  ANGLE  (6,)  LIFT  COEFFICIENT  (PLUNGE)  1C 


PHASE  ANGLE  (02)  LIFT  COEFFICIENT  (PITCH)  1C 


PHASE  ANGLE  (G,)  MOMENT  COEFFICIENT  (PLUNGE) 


CAST  7 


a=0.0° 


Figure  26.  Effect  of  Mach  Number  on  Moment  Coefficient  due  to 
Plunging  for  CAST  7  Airfoil  by  LTRAN2. 


r 


Figure  27.  Effect  of  Mach  Number  on  Moment  Coefficient  due  to 
Pitching  for  CAST  7  Airfoil  by  LTRAN2. 
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decreases.  With  the  increase  in  reduced  frequency,  |C  increases  whereas 
9 2  decreases.  The  slopes  of  all  the  curves  become  much  steeper  when  M  is 
greater  than  approximately  0.7. 

In  Figure  27,  the  magnitude  of  the  moment  coefficient  |C  j  about  the 
quarter  chord  axis  due  to  pitch  about  the  same  axis  increases  as  M  increases. 
The  corresponding  phase  angle  |e^|  first  decreases  and  then  increases  with 
its  minimum  at  M  ;  0.675.  Both  |C  |  and  |e^|  increase  when  the  reduced 
frequency  increases.  The  slopes  for  all  curves  become  much  steeper  when 
M  >  0.7. 

Figures  24  to  27  show  that  the  magnitudes  of  all  the  four  coefficients 
increase  with  the  increase  in  Mach  number.  Such  increase  is  more  rapid  in 
the  higher  Mach  number  region.  The  effect  of  such  changes  will  reflect  on 
the  flutter  behavior. 

4.  Flutter  Results 

Based  on  the  unsteady  aerodynamic  coefficients  obtained  in  Table  3  for 

seven  Mach  numbers,  flutter  analysis  was  performed  for  the  CAST  7  supercritical 

airfoil.  The  emphasis  was  to  investigate  the  effect  of  Mach  number  on  flutter 

speed  for  various  values  of  aeroelastic  parameters. 

Figure  28  shows  the  curves  for  flutter  speed  and  the  corresponding  re¬ 
duced  frequency  versus  Mach  number  for  three  different  values  of  the  position 

of  mass  center  (x  =  0.1,  0.3,  and  0.5,  respectively).  The  values  for  y, 

a 

and  ah  were  assumed  as  100,  0.1,  and  -0.5,  respectively. 

In  this  figure,  the  flutter  speed  increases  as  the  mass  center  approaches 

the  elastic  axis  or  as  x  becomes  smaller.  On  the  other  hand,  the  reduced 

a 

frequency  decreases  as  x&  becomes  smaller. 
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REDUCED  FREQ.(k„)  FLUTTER  SPEED 


All  the  flutter  speed  curves  present  a  dip  phenomenon.  The  dip  appears 
to  be  more  pronounced  and  occurs  at  relatively  lower  Mach  numbers  as  the  mass 
center  approaches  the  elastic  axis. 

Figure  29  shows  the  curves  of  flutter  speed  and  the  corresponding  reduced 

frequency  versus  Mach  number  for  three  values  of  airfoil-air  mass  density 

ratio  (m  =  100,  200,  and  300,  respectively).  The  values  for  the  other  aero- 

elastic  parameters  x  ,  a,  ,  and  w./w  were  assumed  as  0.5,  -0.5,  and  0.1, 

anna 

respectively . 

In  this  figure,  flutter  speeds  are  higher  for  higher  p-values  whereas 
the  corresponding  reduced  frequencies  are  lower  for  higher  p-values.  For 
low  values  of  u,  the  kc-values  corresponding  to  flutter  solution  exceed  0.2 
and  are  beyond  the  limitation  of  LTRAN2.  All  the  three  flutter  curves  pre¬ 
sent  a  dip  phenomenon  with  their  minimum  values  occurring  near  a  Mach  number 
of  0.71. 

Figure  30  shows  the  curves  of  flutter  speed  and  the  corresponding  reduced 

frequency  versus  Mach  number  for  two  values  of  plunge-to-pitch  frequency  ratio 

(ujh/w(i  =  0.1  and  0.2,  respecti vely) .  The  values  for  the  other  aeroelastic 

parameters  y ,  and  a^  were  assumed  as  100,  0.5,  and  -0.5,  respecti vely . 

In  this  figure,  the  flutter  speed  curve  for  =  0.1  is  higher  than 

that  for  uj^/w  =0.2  whereas  the  reduced  frequency  curve  for  the  former  case 

is  lower  than  that  for  the  latter  case.  Flutter  speeds  for  w./w  >  0.2  are 

ha 

not  obtainable  because  the  corresponding  reduced  frequencies  kc  exceed  0.2. 

Both  flutter  curves  present  a  clear  transonic  dip  phenomenon  with  the 
lowest  flutter  speed  occurring  at  M  =  0.71. 

Figure  31  shows  the  curves  of  flutter  speea  and  the  corresponding  re¬ 
duced  frequency  versus  Mach  number  for  two  different  positions  of  elastic 
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REDUCED  FREQ.  (kJ  FLUTTER  SPEED 


REDUCED  FREQ.  (kc)  FLUTTER  SPEED  (U/bw„) 


CAST  7 


CC  s0.0°,  ah  =- 0.5  ,  xa  =  0.5 
=100 


Figure  30.  Effect  of  Mach  Number  on  Flutter  Speed  for  Two  Values  of 
Plunge- to-Pi tch  Frequency  Ratio  for  CAST  7  Airfoil  by 
LTRAN2. 
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REDUCED  FREQ.  (ke)  FLUTTER  SPEED  (U/bw0) 


MACH  NUMBER  (M) 


Figure  31.  Effect  of  Mach  Number  on  Flutter  Speed  for  Two  Positions 
of  Elastic  Axis  for  CAST  7  Airfoil  by  LTRAN2. 
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axis  (ah  =  -0.1  and  -0.5,  respectively).  The  mass  center  was  fixed  at  mid¬ 
chord.  The  other  aeroelastic  parameters,  u  and  were  assumed  as 

100  and  0.1,  respectively. 

In  this  figure,  it  appears  that  the  change  in  flutter  speed  due  to  varia¬ 
tion  of  a^  is  considerably  less  than  the  changes  due  to  variation  of  the  other 
three  aeroelastic  parameters  as  shown  in  Figures  28,  29,  and  30,  respectively. 
The  two  curves  of  corresponding  reduced  frequencies  for  a^  =  -0.5  and  -0.1 
practically  coincide.  Again,  both  flutter  curves  show  a  dip  in  the  neighbor¬ 
hood  of  M  -  0.71 . 

As  explained  in  Reference  29,  the  transonic  dip  is  caused  by  the  compen¬ 
sating  effects  of  the  lift  coefficient  |C£J  due  to  pitch  and  the  position  of 
the  center  of  pressure  ec  (measured  from  the  leading  edge  with  e  being  a  frac¬ 
tion  of  chord  and  c  being  the  full  chord  length).  For  the  CAST  7  airfoil,  plots 
are  presented  in  Figure  32  to  show  such  compensating  effects. 

In  this  figure,  it  is  first  seen  that  |C{J  increases  quite  rapidly  as  M 
becomes  greater  than  approximately  0.68.  Such  increase  has  an  effect  of  re¬ 
ducing  the  flutter  speed.  It  is  then  seen  that  e  remains  as  0.25  (center  of 
pressure  remains  at  1/4-chord)  until  M  '  0.7.  When  M  >  0.7,  e  increases 
rapidly,  i.e.,  the  center  of  pressure  moves  aft.  Such  movement  has  an  effect 
of  increasing  the  flutter  speed.  The  net  effects  of  both  |C  |  and  e  are 
shown  in  the  flutter  curve  with  the  transonic  dip. 
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Figure  32.  Compensate nq  Effects  of  Re  (Cfa)  and  the  Position 
of  Center  of  Pressure  on  Flutter  Speed  for  CAST  7 
Airfoil  by  LTRAN2. 


SECTION  VI 


FLUTTER  ANALYSIS  OF  A  TF-8A  WING  SECTION 


The  TF-8A  supercritical  wing  was  developed  by  NASA  and  it  has  been 
subjected  to  various  aerodynamic  and  aeroelastic  studies  (References  12, 

20,  21,  and  22).  In  this  study,  flutter  analysis  of  the  typical  TF-8A 
wing  section  at  the  65.3%  semispan  station  was  carried  out  by  using  STRANS2 
and  UTRANS2.  Flutter  results  for  the  same  wing  section  were  obtained  by 
Ashley  in  Reference  22  by  using  steady  wind  tunnel  data. 

It  was  noticed  during  flight  tests  of  the  TF-8A  wing  that  because  of 
twist  of  the  wing,  when  corresponding  to  a  zero  mean  angle  of  attack  of  the 
zero-twist  station  at  the  1-g  design  condition,  the  angle  of  attack  at  the 
65.3%  semispan  station  was  -3.0°  (Figure  4  of  Reference  26).  Based  on  the 
wind  tunnel  test  data  from  Reference  26,  a  transonic  dip  was  found  for  the 
TF-8A  wing  section  in  Reference  22.  Therefore,  the  flutter  analysis  for  the 
present  wing  section  was  performed  at  a  -  -3.0°  as  well  as  a  =  0.0°. 

In  this  flutter  analysis,  two  degrees  of  freedom  were  assumed,  plunge 
and  pitch  about  the  1/4-chord  axis.  Seven  Mach  numbers  between  0.7  and  0.8 
were  considered.  Emphasis  was  placed  upon  the  effect  of  Mach  number  on  the 
flutter  speed. 

1.  Airfoil  Configuration 


The  data  for  the  configuration  of  the  TF-8A  wing  section  at  the  65.3% 
semispan  station  were  taken  from  Table  2  of  Reference  26.  Based  on  these 
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measured  data,  the  coefficients  in  Equation  9  for  one-segment  fitting  of  the 
thickness  and  camber  were  provided  by  Dr.  Olsen  of  Flight  Dynamics  Laboratory. 

The  coefficients  for  the  thickness  function  are:  aQ  =  0.3725526; 
a]  =  0.0;  a2  =  -0.5942223;  a3  =  0.9453364;  a4  =  -1.3072624;  and  a&  =  0.5916402. 
The  coefficients  for  the  camber  function  are:  b-j  =  0.0;  b^  =  0.0272615; 
b^  =  -0.2440222;  b^  =  0.6069286;  and  bg  =  -0.3901679.  A  plot  of  the  configura¬ 
tion  is  given  in  Figure  33.  The  maximum  thickness-to-chord  ratio  is  7.9%. 

The  nose  is  relatively  blunt. 

2.  Steady  Pressure  Curves 

The  present  TF-8A  wing  section  has  a  blunter  nose  than  other  supercritical 
airfoils  such  as  MBB  A-3  and  CAST  7.  It  is  desirable  to  evaluate  the  applica¬ 
bility  of  STRANS2  and  UTRANS2  for  this  section. 

Figure  34  shows  two  sets  of  steady  pressure  curves,  one  obtained  by 
STRANS2  and  the  other  by  experiment  in  Reference  26.  The  values  for  Mach 
number  and  mean  angle  of  attack  were  0.5  and  -0.75°,  respectively. 

In  this  figure,  the  two  sets  of  curves  agree  fairly  well  in  overall 
shape.  The  discrepancy  appears  to  be  comparatively  larger  on  the  lower  sur¬ 
face  near  the  trailing  edge.  Similar  discrepancy  was  found  between  LTRAN2  and 
experimental  steady  pressure  curves  obtained  for  the  CAST  7  airfoil  (Figure  22). 

Reasons  for  the  discrepancy  between  the  two  sets  of  curves  can  be  at¬ 
tributed  to:  (1)  viscous  effects  which  are  not  considered  in  STRANS2;  (2) 
small  disturbance  theory  used  in  STRANS2;  and  (3)  possible  losses  of  accuracy 
during  both  the  computational  and  experimental  process. 


75 


AT  65.3%  SEMI  SPAN 


gure  33.  Airfoil  Configuration  for  TF-8A  Wing  Section  at  65.3%  Semispan  Stati 


STEADY  PRESSURE  COEFFICIENT  (c 


Figure  34.  Comparison  between  Steady  Pressure  Curves  by 
STRANS2  and  Experiment  for  TF-8A  Airfoil. 
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Figure  35  shows  the  upper  and  lower  surface  steady  pressure  curves  ob¬ 
tained  for  the  TF-8A  wing  section  for  seven  Mach  numbers  (0.70,  0.74,  0.76, 
0.77,  0.78,  0.79,  and  0.80,  respectively)  at  a  =  -3.0°  by  STRANS2.  The 
steady  pressure  curves  on  the  upper  surface  gradually  increase  with  the  in¬ 
crease  in  Mach  number.  The  lower  surface  steady  pressure  curves  have  shocks 
near  the  leading  edge.  With  the  increase  in  Mach  number,  the  shock  moves 
toward  the  trailing  edge  with  decreasing  strength. 

Figure  36  shows  the  upper  and  lower  surface  steady  pressure  curves  ob¬ 
tained  for  the  TF-8A  wing  section  for  seven  Mach  numbers  (0.70,  0.72,  0.74, 
0.76,  0.78,  0.80,  and  0.82,  respectively)  at  a  =  0.0°  by  STRANS2.  With  the 
increase  in  Mach  number,  the  upper  surface  pressure  curves  change  much  more 
pronounced  than  the  lower  surface  pressure  curves.  For  M  >  0.76,  each  curve 
appears  to  have  two  shocks  on  the  upper  surface,  one  near  the  leading  edge 
and  the  other  near  the  trailing  edge.  The  shock  near  the  leading  edge  de¬ 
creases  in  strength  and  moves  toward  the  trailing  edge  as  Mach  number  in¬ 
creases.  The  shock  near  the  trailing  edge  increases  in  strength  and  moves 
toward  the  trailing  edge  as  the  Mach  number  increases. 

The  shocks  at  the  leading  edge  were  encountered  at  all  Mach  numbers  con¬ 
sidered  here  and  for  both  angles  of  attack  (on  the  lower  surface  at  a  =  -3.0° 
and  c n  the  upper  surface  at  u  =0°).  They  seem  to  be  due  to  the  effect 
of  the  blunt  nose  of  the  airfoil. 

3.  Unsteady  Aerodynamic  Coefficients 

Four  unsteady  aerodynamic  coefficients,  C„r,  C.  ,  C  r,  and  C  ,were  ob- 

x-6  y.a  mo  ma 

tained  for  the  TF-8A  wing  section  by  UTRANS2.  The  airfoil  was  pitched  about 
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Figure  35.  Effect  of  Mach  Number  on  Steady  Pressure  Curves  for  TF-8A 
Airfoil  at  a  =  -3.0°  by  STRANS2. 
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Figure  36.  Effect  of  Mach  Number  on  Steady  Pressure  Curves  for  TF-8A 
Airfoil  at  a  =  0.0°  by  STRANS2. 
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the  quarter  chord  at  various  Mach  numbers  for  two  angles  of  attack 


{  «  --  -3.0"  and  «  =  0°). 

Mach  numbers  considered  were  0.7,  0.74,  0.76,  0.77,  0.78,  0.79,  and 
0.80  for  u  =  -3.0°  and  0.7,  0.72,  0.74,  0.76,  and  0.78  for  a  =  0.0°. 

The  reduced  frequencies  assumed  for  both  cases  were  0.0,  0.05,  0.10, 

0.15,  and  0.20.  Table  4  and  Table  5  give  the  four  unsteady  coefficients 
obtained  at  a  =  -3.0°  and  a  =  0.0°,  respectively.  These  four  coefficients 
and  ti.  ■  corresponding  phase  angles  are  also  plotted  against  Mach  number  in 
Figures  37,  38,  39,  and  40,  respectively. 

In  Figure  37,  the  magnitude  of  the  lift  coefficient  due  to  plunge  jC„rJ 
for  a  =  -3.0°  increases  very  slightly  as  the  Mach  number  increases.  Tne 
curve  shows  a  bump  in  the  neighborhood  of  M  *  0.78.  The  corresponding  phase 
angle  decreases  slightly  as  M  increases.  As  the  reduced  frequency  in¬ 
creases,  |C  j  increases  whereas  6^  decreases.  Also  shown  in  the  same  figure 
are  those  curves  for  a  =  0.0°.  In  this  case,  slightly  flatter  bumps  in  the 
;C,6|  curves  appear  in  the  neighborhood  of  M  =  0.76. 

In  Figure  38,  the  magnitude  of  the  lift  coefficient  jCfJ  due  to  pitcn 

about  the  quarter  chord  axis  for  a  =-3.0°  increases  as  M  increases.  Again, 

bumps  are  present  in  the  neighborhood  of  M  =  0.78.  The  corresponding  phase 
angle  1621  increases  as  M  increases.  As  the  reduced  frequency  increases, 
jCj  j  decreases  whereas  je^j  increases.  Also  shown  in  the  same  figure  are 
those  curves  for  u  =  0.0°.  The  bumps  are  much  milder  and  appear  at  lower  M 
than  those  for  a  =  -3.0°. 

In  Figure  39,  the  magnitude  of  the  moment  coefficient  |C  ,1  about  the 

quarter  chord  axis  due  to  plunge  for  a  =  -3.0°  increases  as  M  increases. 

The  corresponding  phase  angle  decreases  as  M  increases.  As  the  reduced 
frequency  increases,  IC^I  increases  whereas  6^  does  not  show  any  clear 
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Table  4. 


Aerodynamic  Coefficients  for  TF-8A  Airfoil  for 
Various  Mach  Numbers  at  a  =  -3.0°  by  UTRANS2 
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Table  5.  Aerodynamic  Coefficients  for  TF-8A  Airfoil  for 
Various  Mach  Numbers  at  a  =  0.0°  by  UTRANS2 
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Figure  37.  Effect  of  Mach  Number  on  Lift  Coefficient  due  to 
Plunging  for  TF-8A  Airfoil  by  UTRANS2. 
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trend.  Also  shown  in  the  same  figure  are  those  curves  for  a  =  0.0°. 

These  curves  show  a  similar  trend  as  those  for  a  =  -3.0°  but  are  smoother. 

In  Figure  40,  the  magnitude  of  the  moment  coefficient  |Cm  |  about  the 
quarter  chord  axis  due  to  pitch  about  the  same  axis  for  a  =  -3.0°  increases 
as  M  increases.  The  corresponding  phase  angle  decreases  as  M  increases. 

As  kc  increases,  |CmJ  increases  when  M  <  0.785  approximately  and  decreases 
when  M  >  0.785  approximately.  As  kc  increases,  9 ^  increases  when  M  >  0.72. 
Also  shown  in  the  same  figure  are  those  curves  for  a  =  0.0°.  Both  sets  of 
curves  show  a  similar  trend.  However,  the  crossover  of  the  |C  |-curves  for 
a  =  0.0° takes  place  at  M  1  0.75  and  the  crossover  of  the  e^-curves  takes 
place  at  M  =  0.78. 

4.  Flutter  Results 

Based  on  the  unsteady  aerodynamic  coefficients  obtained  for  various 
Mach  numbers  in  Tables  4  ai.d  5  for  a  =  -3.0°  and  0.0°,  respectively,  flutter 
analysis  was  carried  out  for  the  TF-8A  wing  section.  The  emphasis  was  to 
investigate  the  effect  of  Mach  number  on  flutter  speed  for  various  values  of 
aeroelastic  parameters. 

Figure  41  shows  curves  of  flutter  speed  and  corresponding  reduced 

frequency  versus  Mach  number  at  a  =  -3.0°  for  four  values  of  the  position  of 

mass  center  (x  =  0.3,  0.4,  0.5,  and  0.6).  The  values  for  y,  oi./w  ,  and  a. 

a  nan 

were  assumed  as  100,  0.1,  and  -0.5,  respectively. 

In  this  figure,  the  flutter  speed  curves  shift  to  higher  levels  as  the 
mass  center  moves  toward  the  elastic  axis  or  as  xq  decreases.  The  correspond¬ 
ing  reduced  frequency  curves  show  a  reversed  trend. 
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Figure  41.  Effect  of  Mach  Number  on  Flutter  Speed  for  Four 
Positions  of  Mass  Center  for  TF-8A  Airfoil  by 
STRANS2/UTRANS2. 


For  a  =  -3.0°,  all  the  flutter  curves  show  a  dip  in  the  neighborhood 
of  M  =  0.78.  The  dip  appears  to  be  more  pronounced  when  the  mass  center  is 
closer  to  the  elastic  axis. 

Also  shown  in  Figure  41  are  the  flutter  results  for  a  =  0.0°  for  com¬ 
parison.  No  obvious  dips  are  present  in  these  flutter  curves.  However, 
it  appears  that  the  flutter  speeds  are  lowest  in  the  neighborhood  of  M  =  0.72. 

Figure  42  shows  the  flutter  results  for  both  cases:  a  =  -30°  and 
a  -  0.0°.  Three  y-values  of  100,  200,  and  300,  respectively,  were  considered. 
The  values  selected  for  xa>  a^,  and  ^/wa  were  0.5,  -0.5,  and  0.1,  respectively. 

In  this  figure,  the  flutter  speed  curves  are  higher  at  higher  altitude  or 
at  higher  y-values.  The  corresponding  reduced  frequency  curves  show  a  re¬ 
versed  trend. 

All  three  flutter  curves  for  a  =  -3.0°  present  a  dip  in  the  neighborhood 
of  M  =  0.78°.  The  three  curves  for  a  =  0.0°  do  not  show  obvious  dips.  These 
three  curves  are  wavy  and  each  has  two  lower  points  in  the  neighborhood  of 
M  =  0.72  and  0.76,  respectively.  However,  the  curves  obviously  show 
increasing  trend  when  the  Mach  number  goes  beyond  M  =  0.76. 

Figure  43  shows  the  flutter  results  for  both  a  =  -3.0°  and  0.0°  and  for 
=  0.1  and  0.2,  respectively.  The  values  selected  for  x^,  a^,  and  y  were 
0.5,  -0.5,  and  100,  respectively. 

In  this  figure,  the  flutter  speed  curves  for  u^/uj^  =  0.1  are  higher  than 
those  for  wh/u>a  =  0.2.  The  corresponding  reduced  frequency  curves  have  an 
opposite  relationship  with  the  two  frequency  ratios. 

Again,  the  two  flutter  curves  for  a  =  -3.0°  present  a  dip  in  the  neighbor¬ 
hood  of  M  =  0.78.  The  two  flutter  curves  for  a  =  0.0°  do  not  show  obvious 
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dips  except  two  very  mild  lower  waves,  one  near  M  =  0.72  and  the  other  near 
M  =  0.76.  However,  the  flutter  speeds  do  tend  to  rise  obviously  when  the 
Mach  number  goes  beyond  M  =  0.76. 

Figure  44  shows  the  flutter  results  for  both  a  =  -3.0°  and  0.0°  and  for 

ah  =  -0.1  and  -0.5,  respectively.  The  mass  center  was  fixed  at  mid-chord  as 

the  elastic  axis  moved.  The  values  for  u  and  w./w  were  assumed  as  100  and 

h  a 

0.1 ,  respectively. 

For  either  case  (a  =  -3.0°  or  0.0°),  the  flutter  speed  curves  for 
a^  =  -0.1  and  a^  =  -0.5  are  not  much  different.  Such  differences  are  even  less 
obvious  between  the  curves  for  the  corresponding  reduced  frequencies. 

As  also  shown  in  the  previous  two  figures,  the  two  flutter  speed  curves 
for  a  =  -3.0°  exhibit  a  transonic  dip  in  the  neighborhood  of  M  =  0.78.  The 
two  flutter  speed  curves  for  a  =  0.0°  do  not  exhibit  obvious  dips.  However, 
the  two  curves  do  show  obvious  rising  tendency  when  M  becomes  greater  than  0.76. 

As  discussed  in  Reference  29,  the  transonic  dip  phenomenon  may  be  physi¬ 
cally  explained  by  looking  at  the  compensating  effects  of  the  increase  in  mag¬ 
nitude  of  the  lift  coefficient  j |  due  to  pitch  which  tends  to  decrease  the 
flutter  speed  and  the  aft-movement  of  the  center  of  pressure  which  tends  to  in¬ 
crease  the  flutter  speed.  When  the  airfoil  goes  supersonic,  the  center  of 
pressure  moves  to  the  mid-chord  on  the  airfoil  and  flutter  speed  becomes  much 
higher. 

The  real  part  of  lift  coefficient  (Cfa),  the  position  nf  center  of 
pressure  e,  and  the  flutter  speed  are  plotted  against  Mach  number  in  Figure 
45  for  the  TF-8A  section.  The  position  of  the  center  of  pressure  is  defined 
as  at  a  distance  e  times  c  aft  of  the  leading  edge. 
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For  the  case  a  =  -3.0°,  the  lift  coefficient  curve  shows  an  obvious 
peak  at  M  =  0.78  which  strongly  suggests  a  dip  in  flutter  speed.  The 
center  of  pressure  plot  shows  a  gradual  and  steady  aft-shift  from  1/4-chord 
axis  when  M  becomes  higher  than  0.76.  Such  movement  suggests  a  comparable 
increase  in  flutter  speed.  The  flutter  speed  curve  is  a  result  of  these  two 
effects. 

For  a  =  0.0°,  the  lift  coefficient  curve  suggests  decrease  in  flutter 
speed  with  no  obvious  dip  whereas  the  e-curve  suggests  an  increase  in  flut¬ 
ter  speed  earlier  and  sharper  than  that  for  a  =  -3.0°. 
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SECTION  VII 


CONCLUDING  REMARKS 


The  purpose  of  this  research  was  to  study  the  flutter  behavior  of 
three  supercritical  airfoils:  MBB  A-3,  CAST  7,  and  TF-8A,  through  the  use 
of  LTRAN2  and  STRANS2/UTRANS2  so  that  some  indications  or  conclusions  may 
be  drawn  in  the  aspects  such  as:  (1)  applicability  and  limitation  of  the 
two  codes;  (2)  trends  and  flutter  behavior  of  the  subject  airfoils  for 
the  parameters  considered;  and  (3)  correlation  and  interpretation  of  the 
results  for  various  cases.  In  the  interest  of  simplicity,  the  airfoil 
geometries  were  represented  by  the  simplest  possible  one-segment  fits  with 
a  minimum  of  constraints.  More  elaborate  descriptions  may  be  necessary 
for  highly  cambered  airfoils. 

For  the  MBB  A-3  airfoil,  both  codes  were  used  to  investigate  the  effect 
of  small  mean  angles  of  attack  on  its  flutter  behavior  at  the  design  Mach 
number  of  0.765. 

For  the  CAST  7  airfoil,  LTRAN2  was  used  to  investigate  the  effect  of  Mach 
number  on  its  flutter  behavior  at  zero  mean  angle  of  attack. 

For  the  TF-8A  wing  section  at  the  65.3%  semispan  station,  STRANS2/UTRANS2 
were  used  to  investigate  the  effect  of  Mach  number  on  its  flutter  behavior 
at  a  =  0.0°  and  -3.0°  (In  Reference  26,  when  the  TF-8A  wing  was  tested  in  the 
1-g  design  cruise  condition  with  zero  angle  of  attack  at  the  zero-twist  sta¬ 
tion,  the  65. 3%-section  had  an  a-value  of  -3.0°). 

As  a  result  of  this  study,  the  following  concluding  remarks  are  made. 
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a.  In  general,  flutter  results  obtained  by  both  codes  for  the  MBB 
A- 3  airfoil  at  M  =  0.765  appear  to  show  a  similar  trend.  For  the  para¬ 
meters  assumed  (u  =  100,  200,  300;  x^  =  0.1,  0.3,  0.5;  a^  =  -0.1,  -0.5; 
and  wh/wa  =  0.1),  most  of  the  flutter  speed  plots  show  a  trend  that  in¬ 
creases  with  angle  of  attack.  This  is,  however,  not  the  case  for  those 

obtained  by  LTRAN2  at  w./w  =  0.2. 

n  a 

b.  For  the  same  flow  conditions,  LTRAN2  gives  stronger  shocks  than 
STRANS2. 

c.  Because  of  the  limitations  in  both  LTRAN2  and  STRANS2/UTRANS2 , 

0.8°  and  1.3°  are  approximately  the  highest  a  the  two  respective  programs 
can  treat  for  the  MBB  A- 3  cases  considered. 

d.  For  the  CAST  7  airfoil,  the  steady  pressure  curves  obtained  by 
LTRAN2  at  M  =  0.602  and  a  =  0.5°  compare  fairly  well  with  those  obtained 
from  wind  tunnel  test  (Reference  11)  in  overall  shape.  However,  the  LTRAN2 
coefficients  are  higher  than  those  from  test. 

e.  For  all  the  parameters  assumed  for  the  CAST  7  airfoil,  the  various 
plots  of  the  flutter  speed  curves  show  a  dip  in  the  neighborhood  of  M  =  0.705. 

f.  For  the  TF-8A  wing  section,  the  steady  pressure  curves  obtained  by 
STRANS2  at  M  =  0.5  and  a  =  0.75°  compare  fairly  well  with  experimental  data 
(Reference  26). 

g.  For  all  the  parameters  assumed  for  the  TF-8A  section,  the  various 
plots  of  the  flutter  speed  curves  show  a  dip  in  the  neighborhood  of  M  =  0.78 
when  a  =  -3.0°.  The  plots  for  the  case  o  =  0.0°  do  not,  however,  show  ob¬ 
vious  dips  other  than  a  drop  in  flutter  speed  over  the  range  0.71  £  M  < 


h.  The  reasons  for  the  discrepancies  between  the  steady  pressure 
curves  obtained  by  computations  and  experiments  in  Figures  22  and  34  may 
be  attributed  to  the  reasons  that  (a)  viscous  effects  have  not  been  in¬ 
cluded  in  both  computer  codes,  (b)  small  disturbance  theory  is  used  in 
both  codes,  and  (c)  inaccuracies  may  have  been  introduced  in  both  the  nu¬ 
merical  and  experimental  procedures. 

i.  Numerical  convergence  problems  are  encountered  wher.  the  shock 
grows  with  Mach  number.  For  the  cases  studied,  upper  limits  of  Mach  number 
have  been  found  by  both  codes. 

j.  Due  to  the  low-frequency  approximation  used  in  LTRAN2,  the  reduced 

frequencies  k£  are  limited  to  be  not  higher  than  0.2.  Although  the  low- 

frequency  approximation  is  not  used  in  UTRANS2,  convergence  difficulty  is 

usually  encountered  when  k  >  0.2. 

c 

k.  For  a  more  complete  understanding  of  the  transonic  flutter  be¬ 
havior  of  supercritical  airfoils,  the  present  computer  codes  have  to  be  im¬ 
proved  to  be  applicable  to  a  wider  range  of  flow  and  aeroelastic  parameters. 
The  limitations  on  Mach  number,  angle  of  attack,  inviscid  flow,  small  dis¬ 
turbances,  and  reduced  frequency,  etc.,  must  be  relaxed.  Such  developments 
have  been  under  way  (References  18  and  19). 

l.  More  developments  in  3-D  transonic  codes  and  their  applications  to 
flutter  analysis  of  full  wings  are  needed.  Some  preliminary  3-D  flutter 
analysis  using  TDSTRN  and  TDUTRN  was  given  in  Reference  30. 

m.  The  present  results  may  be  of  use  to  those  who  are  interested  in 
the  transonic  codes  and  flutter  behavior  of  supercritical  airfoils.  They 
may  also  serve  as  a  comparative  basis  for  experimental  results.  On  the  other 
hand,  more  experimental  unsteady  pressure  and  flutter  data  are  needed. 
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